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We studied the energy spectrum and the large-scale fluctuation of the X-ray background with the ASCA 
' GIS instrument based on the ASCA Medium Sensitivity Survey and Large Sky Survey observations. A 

total of 91 fields with Galactic latitude |6| > 10° were selected with a sky coverage of 50 deg 2 and 4.2 
Ms of exposure. For each field, non X-ray events were carefully subtracted and sources brighter than 
' ~2x 10~ 13 erg cm -2 s" 1 (2-10 keV) were eliminated. Spectral fits with a single power-law model for 

the individual 0.7-10 keV spectra showed a significant excess below ~ 2 keV, which could be expressed 
by an additional thermal model with kT ~ 0.4 keV or a steep power-law model with a photon index of 
psoft ^ g rp^g q g_2 k e y intensities of the soft thermal component varied significantly from field to field 
■ by 1 a = 52TJ %, and showed a maximum toward the Galactic Center. This component is considered 

| to be entirely Galactic. As for the hard power-law component, an average photon index of 91 fields 

£N| i was obtained to be r hard = 1.412 ± 0.007 ± 0.025 and the average 2-10 keV intensity was calculated as 

^hard _ (6.38 ±0.04 ±0.64) x 10" 8 erg cm -2 s" 1 sr" 1 (1 a statistical and systematic errors). The Galactic 
component is marginally detected in the hard band. The 2-10 keV intensities show a 1 <r deviation of 
6.49^061 %, while deviation due to the reproducibility of the particle background is 3.2%. The observed 
deviation can be explained by the Poisson noise of the source count in the f.o.v. (~ 0.5 deg 2 ), even assuming 
i-^ ' a single log AT- log<S relation on the whole sky. Based on the observed fluctuation and the absolute intensity, 

an acceptable region of the log A- logS 1 relation was derived, showing a consistent feature with the recent 
Chandra and XMM-Newton results. The fluctuation of the spectral index was also examined; it implied a 
\ large amount of hard sources and a substantial variation in the intrinsic source spectra (Tg ~ 1.1 ± 1.0). 

c/3 . 

, 1. Introduction with a uniform distribution of sources in Euclidean space, 

namely, A(> S) oc S" 15 (Piccinotti et al. 1982; Hayashida 
. £h ! The Cosmic X-ray Background, hereafter called CXB, et al. 1992; Ueda et al. 1999a, b). At the present time, 
dominates the hard X-ray surface brightness (Fabian, the 2-10 keV CXB has been almost resolved into discrete 
Barcons 1992), and was discovered with the Aerobee sources, leaving at most 10-20% at the faintest flux limit, 
sounding rocket carrying three Geiger counters (Giacconi The main uncertainty of the resolved fraction lies in the 
et al. 1962) along with the first extra-solar X-ray source, measurement of the total intensity of the CXB, which re- 
Sco X-l. The effort of resolving the CXB into discrete quires a well-calibrated detector with a low internal back- 
sources has produced a log A- logS* relation, which gives ground, as well as large sky coverage, 
the number density (A) of discrete sources above a cer- It also became clear that the emission consists of 
tain flux level (S). This study made a major advance re- roughly two components. The hard-band emission in ~ 2- 
cently with the advent of the Chandra X-ray Observatory. 10 keV shows almost an isotropic distribution in the sky, 
The detection limit of Chandra is down to S ~ 10" 16 and the spectrum is characterized by a power-law spec- 
erg cm" 2 s" 1 in the 2-10 keV band, and some of the trum with a photon index of T ~ 1.4 (e.g. Marshall et 
sources are too faint for optical identifications to be made. al. 1980: 3-15 keV; Gendreau et al. 1995: 1-10 keV; 
In this energy range, more than half (60-90%) of the Miyaji et al. 1998: 1-10 keV; Vecchi et al. 1999: 1- 
CXB emission has been resolved into discrete sources 8 keV). The survey observations by Uhuru (2-7 keV) 
(e.g. Mushotzky et al. 2000 for the SSA13 field; Brandt and HE AO 1 A-2 (~ 2-10 keV) showed that the dipole 
et al. 2001 for the Chandra Deep Field North; Tozzi et amplitude after removing the Galactic contribution were 
al. 2001 for the Chandra Deep Field South). Moreover, 0.61 ±0.26% (Protheroe et al. 1980) and - 0.5% (Shafer, 
recent Chandra and XMM-Newton results indicate a turn- Fabian 1983), respectively. Below ~ 2 keV, the spectrum 
over of around S ~ 2 x 10~ 14 erg cm -2 s" 1 (Campana becomes steeper with T~ 2.1 (e.g. Hasinger 1992 at ~ 1 
et al. 2001; Hasinger et al. 2001; Baldi et al. 2002). keV). The R OSA Tall-sky survey observation showed sig- 
Above this flux level, the relation seems to be consistent nificant structures in the soft X-ray background which 



2 



A. Kushino et al. 



are correlated mainly with the Galactic distribution of 
the hot gas (Snowden et al. 1995). As for the absolute 
intensity of the CXB, there is fairly large uncertainty 
among the measurements; namely, the reported intensi- 
ties at 1 keV are 13.4 ±0.3 (Hasinger 1992), 9.4 ±0.4 
(Gendreau et al. 1995) and 10.4+^ (Parmar et al. 1999) 
photon kcV" 1 cm -2 s _1 sr _1 . Barcons ct al. (2000) com- 
piled previous measurements of the CXB intensity with 
the ASCA, BeppoSAX and ROSAT and pointed out that 
these differences are primarily caused by the cosmic vari- 
ance, i.e. spatial variation of source count due to the lim- 
ited solid angle of sky coverage, and in some part resulted 
from instrumental cross-calibration errors and subtraction 
process of the Galactic contribution. 

X-ray surveys with optical follow-up observations (e.g. 
Schmidt ct al. 1998; Akiyama et al. 2000; Ishisaki et 
al. 2001) have identified many active galactic nuclei 
(AGNs) and provided information about their spectral 
evolution. There is a significant discrepancy between the 
spectral indices of the CXB above ~ 2 keV (r ~ 1.4) 
and of type-1 Seyfert galaxies (r ~ 1.7), which led to a 
model that a significant fraction of AGN is heavily ab- 
sorbed (Madau et al. 1994; Comastri et al. 1995; Gilli 
et al. 2001). In fact, for the sources with S > 10 -15 
erg cm~ 2 s _1 , the XMM-Newton satellite revealed that 
around half of them are strongly absorbed in their spectra 
(Hasinger et al. 2001). In the flux range S ~ lO^-lO" 11 
erg cm -2 s _1 , a spectral survey from ASCA showed that 
the average spectral slope becomes progressively harder 
as the sources become fainter, and the observed photon 
index of 2.1 at around 10 -11 erg cm -2 s _1 turns to 1.6 
at around 10 -13 erg cm -2 s" 1 in the 0.7-10 keV band 
(Ueda et al. 1999b). A precise determination of the spec- 
tral shape of the CXB and an investigation of its field-to- 
field difference would provide rich information to under- 
stand the spectral evolution of AGNs taking place around 
the flux of 10 -13 erg cm -2 s _1 . 

The ASCA Medium-Sensitivity Survey (AMSS) is a 
serendipitous survey, consisting of 368 fields which can 
be regarded as a random sampling of the CXB. Ueda et 
al. (2001) produced a catalog of 1343 X-ray sources based 
on the AMSS, and gave a fairly tight log AT- logS relation 
in the flux range above 7 x 10 -14 erg cm~ 2 s _1 (Ueda et 
al. 1999b). Also, combined with the ASCA Large Sky 
Survey (Ueda ct al. 1998, 1999a), the main aim of the 
present study is to determine the absolute intensity and 
the spectral shape of the CXB as precisely as possible. 
The data after the point-source elimination formed the 
basis for the present study. It would also enable us to 
constrain the log A- \ogS relation or the intrinsic spec- 
tral distribution of the constituent sources through the 
fluctuation analysis, as well as to look into the large-scale 
distribution of the CXB. The GIS system (Ohashi et 
al. 1996; Makishima et al. 1996) consists of two gas imag- 
ing spectrometers, which have well-studied and low inter- 
nal background, simultaneously providing a large (^ 40' in 
diameter) field of view (f.o.v.), when compared with the 
XMM-Newton PN and MOS CCDs and Chandra ACIS 
CCDs. These properties are powerful for studying the 



spectral features of the CXB. 

The remainder of this paper is organized as follows. In 
section 2, we discuss the selection criteria of the fields 
employed in the present paper for the CXB analysis and 
explain how energy spectra are extracted. In section 3, 
we present the fitting results. In section 4, the large-scale 
distributions in the spectral parameters are investigated. 
In section 5, using simulations which take into account 
the complicated response of ASCA, we evaluate the CXB 
intensity and constrain the log A- log S relation and the 
spectral distribution of sources from the observed fluctu- 
ation. In section 6, we discuss implications of the results, 
and the summary will be given in section 7. 

2. Observation and Analysis 

The present study was carried out with a GIS instru- 
ment consisting of GIS 2 and GIS 3 sensors (Ohashi et 
al. 1996; Makishima et al. 1996). The detectors had a wide 
(~ 40' in diameter) f.o.v. and a moderate spatial resolution 
of ~ 3' (Half Power Diameter) combined with the X-Ray 
Telescope (XRT; Serlemitsos et al. 1995). Our primary 
goal of the analysis was to determine the parameters of the 
CXB for each direction which cannot be resolved into dis- 
crete sources with the GIS, namely, the unresolved CXB. 
Due to the low surface brightness of the unresolved CXB, 
we had to take into account several issues that would not 
matter very much in point-source studies. First of all, the 
sample fields must be carefully selected in order to obtain 
reasonable constraints on the parameters. The second one 
is to minimize and precisely estimate the amount of fake 
signals caused mainly by charged particles, which we call 
the Non X-ray Background (NXB). The third one is to 
eliminate resolved sources in the GIS f.o.v. The last one 
is the instrumental response for the unresolved CXB. 

2.1. Selection of the Fields 

Our sample fields were primarily selected from the first 
AMSS fields (Ueda et al. 2001), which consisted of 368 
combined fields observed between 1993 May and 1996 
December with the Galactic latitude \b\ > 10°. We also in- 
cluded the field of the ASCA Large Sky Survey (LSS; Ueda 
et al. 1998, 1999a), dividing the LSS field into four, label- 
ing them a-LSS, b-LSS, c-LSS, and d-LSS, from south- 
west to north-east. 

In our prompt sample, several fields were composed of 
multiple pointing observations, covering a larger sky area 
than that of a single pointing. Furthermore, most of the 
AMSS fields were pointed to certain targets and often 
contained rather bright X-ray sources to study the CXB. 
Since the Point Spread Function (PSF) of the XRT has 
long outskirt, approximated by 

PSF(r) = exp (-^) + 0.0094 exp (-^) , (1) 

a bright target source would contaminate the whole de- 
tector area (r < 20'). Wc therefore selected fields based 
on the following criteria: 

1. Total exposure time summed up over all pointings 
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is longer than ~ 20 ks. 

2. No known extended X-ray sources is present in the 
field, e.g., nearby galaxies, star- forming regions, or 
supernova remnants. 

3. The mean count rate including the NXB is lower 
than 0.3 count s _1 per one GIS sensor in 0.7-10 
keV, accumulated within a radius of 22 mm (~ 22') 
from the optical axis of each sensor. For multi- 
pointing observations, we accept fields when at least 
one pointing fulfills this condition. 

4. The remained sky area after the source elimination 
("Area 2" column of table 1) is more than two thirds 
of the original sky area ( "Area 1" ) . 

Exactly 100 fields satisfied these conditions. 

Since the XRT allows some fraction of photons from 
outside of the f.o.v. to pass through (stray light; Tsusaka 
et al. 1995), it is a significant problem in the CXB study 
with ASCA. Figure 1 represents the 0.7-10 keV count rate 
of the Crab nebula pointed at various offset angles from 
the optical axis of each GIS sensor. When the Crab nebula 
was placed at ~ 1° away from the optical axis, namely far 
outside of the GIS f.o.v., its count rate was still larger than 
10 count s -1 sensor -1 in 0.7-10 keV. Since the CXB gave 
only ~ 0.12 count s -1 sensor -1 in the same energy band, 
the stray flux needed to be carefully examined. In order 
to avoid any influence of the stray light in the present 
analysis, we further selected those fields which had no 
bright sources around the f.o.v. We mainly consulted with 
the flOS^TAll-Sky Survey Bright Source Catalog (RASS- 
BSC; Voges et al. 1999), which lists sources with a flux 
limit of 2.8 x 10~ 13 erg cm -2 s -1 , assuming a power-law 
photon index of T = 2 in the 0.1-2.4 keV energy band. 
However, the sensitive energy range of RO SAT is limited 
to below ~ 2 keV and sources with hard spectra or strong 
absorption would be missed. The HEAO 1 A-l X-ray 
source catalog (Wood et al. 1984) was also used to follow 
up such cases. HEAO 1 A-l is sensitive above ~ 2 keV, 
and the flux limit is 4.78 x 10 -12 erg cm -2 s -1 in 2-10 
keV for a Crab-like spectrum. 

Based on these catalogs, the distance distribution of 
the cataloged sources from the centers of the selected 100 
fields was determined as shown in figure 2. In this figure, 
the source intensities are plotted as a function of the off- 
set angle. The solid curve indicates a level where sources 
above it could contribute to the CXB intensity accumu- 
lated within the whole GIS f.o.v. (r < 20 mm ~ 20') by 
more than 2.5%. This estimation is based on a calibration 
observation of the Crab nebula at various offset angles, 
which is indicated in figure 1 by the solid line. Using the 
ROSAT sources, we excluded 9 fields in which there were 
sources above the 2.5% line outside the GIS f.o.v. (r > 20 
mm ~ 20'). We did not exclude fields if bright sources are 
inside the f.o.v., i.e. in the hatched region in figure 2, since 
they would be eliminated in the source-elimination pro- 
cess (see subsection 2.3). As for the HEAO 1 sources, any 
fields were excluded if sources brighter than 1.0 x 10 -10 
erg cm -2 s _1 were present within r < 4°. This is because 
the positional errors of the HE A O 1 sources are quite large 
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(~ 1°). One field was picked up by the HEAO -/-source 
selection; however, it had already been marked by the 
ROSA T-source selection. We therefore excluded 9 fields 
in total. One of the LSS fields (a-LSS) was rejected here. 

After these field selections, we ended up with 91 fields 
whose sky distribution is shown in figure 3 together with 
the original AMSS and LSS fields. The names and coor- 
dinates of the fields are listed in table 1 . Among them, 33 
fields consist of multi-pointing observations. The exposure 
time of the 91 (88) fields amounts to 4240.54 (3844.27) 
ks and the covering sky area extends 49.53 (42.77) deg 2 , 
where the numbers in parentheses denote those with- 
out the three LSS fields. According to the AMSS and 
the LSS catalogs, the brightest source in the 91 fields is 
1AXGJ 085449+2006 with 0.7-7 keV flux of 5.89 x 10 -12 
erg cm" 2 s" 1 in OJ 287 field, or 1AXG J004847+3157 
with 2-10 keV flux of 7.58 x 10 -12 erg cm -2 s" 1 in Mrk 348 
field, depending on the chosen energy band. Both sources 
correspond to the target, itself, i.e. OJ 287 or Mrk 348. 
Among the 88 AMSS fields, there are 184 (141) sources 
with their 0.7-7 (2-10) keV flux brighter than 2 x 10 -13 
erg cm -2 s -1 in each energy band in the AMSS catalog. 

2.2. Data Screening and NXB Subtraction 

We screened all of the GIS events detected in each sam- 
ple field, first employing the standard event selection cri- 
teria: (1) the GIS should be in the nominal observation 
mode, i.e. the PH normal mode with the nominal bit as- 
signment, the spread discriminator turned on, and the 
satellite not in the South Atlantic Anomaly (SAA) where 
the high-voltage supply of the GIS is switched off; (2) the 
elevation of f.o.v. should be > 5° (or 25°) above the night 
(or sunlit) Earth rim; (3) the geomagnetic cut-off rigidity 
(COR) should be > 6 GV. We performed further event 
screening utilizing the GIS monitor counts (HO + H2; see 
Ohashi et al. 1996; Makishima et al. 1996), Radiation Belt 
Monitor (RBM) count, and COR. This additional selec- 
tion improved the reproducibility of the NXB estimation, 
mainly by rejecting data affected by sporadic increases 
in the NXB count, which were due probably to the con- 
centration of charged particles on the satellite orbit. A 
stricter condition for the RBM count was applied around 
the SAA and Hawaii. Details of the screening procedure 
are described in Ishisaki (1996) and Ishisaki et al. (1997). 

The obtained event list contains not only the CXB pho- 
tons, but also the NXB events, which must be carefully 
subtracted. The corresponding spectra or images of the 
NXB were separately created, using observations when 
the XRT was pointed to the night Earth. Under the 
same condition as described in the previous paragraph, 
we processed seven years of night-Earth data from 1993 
June to 2000 May, with a net exposure time of 4880.43 
ks. We further sorted both the observed data and the 
night-Earth data into six intervals of the HO + H2 moni- 
tor count in 5 count s _1 steps, and performed NXB sub- 
traction for each interval. This is because the HO + H2 
monitor count, which is tightly correlated with the NXB 
count (Makishima et al. 1996), can be used as a good NXB 
indicator during on-source observations. 
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There is also a long-term change of the NXB count rate, 
which is shown in figure 4. The count rate gradually in- 
creased for the first four years, and then turned to decreas- 
ing, peaking at around 1997-1998. This long-term change 
was presumably caused by a gradual drop in the satellite 
altitude and the cycle of the solar activities with a period 
of ~ 11 yrs. During the quiescent state of the solar ac- 
tivities, which is at the minimum around 1996-1997, the 
atmosphere shrunk so that the density of charged particles 
increased. Since the observation date of our sample fields 
ranged from 1993 June to 1996 December, there could be 
about a ±10% error in the NXB estimation unless the 
long-term trend were corrected. We fitted the trend by a 
fourth-order polynomial, which was used for the correc- 
tion. 

We estimated that the NXB had been reproduced by 
an accuracy of 3% using the HO + H2 monitor count and 
corrected for the long-term trend. Therefore, the system- 
atic uncertainty in the NXB subtraction affects the CXB 
intensity by at most 10%, even at ~ 7 keV where the NXB 
dominates the CXB (see figure 7). 

2.3. Source Elimination 

The CXB intensity can be defined by equation (6) in 
subsection 5.1, even including the contribution from very 
bright sources, namely, by setting S ^ oo in equation 
(6) in subsection 5.1. On the other hand, because ac- 
tual observations are usually biased to the faint so-called 
blank sky, the observed intensity needs a certain correc- 
tion. This correction depends on the observed flux range, 
S < So, where it is assured that only the sources with their 
fluxes fainter than So exist in a sample field. Therefore, 
it is essentially required to eliminate resolved sources on 
the source flux basis. We employed a certain fixed value 
of S ~ 2.0 x 10~ 13 erg cuT 2 s" 1 (2-10 keV), and defined 
that sources brighter than Sq were resolved and that the 
remained X-ray emission was the unresolved CXB. This 
procedure is also important in order to overcome the se- 
lection bias of the AMSS fields, because they are origi- 
nally intended to observe specific targets. Since the AMSS 
and the LSS catalogs are constructed on the source sig- 
nificance basis, some sources with low significances but 
possible high flux or extended emission could be missed. 
We therefore developed an original method for eliminating 
sources by slightly modifying the source detection proce- 
dure of the AMSS. Throughout this paper, we call this 
method as "source elimination" , which is described in the 
following paragraphs. 

The resolved sources were eliminated using the observed 
data, itself, in each field. We first constructed a mosaic 
flat-field image M FLAT in the sky coordinates. This is a 
probability map that is proportional to the intensity of 
resolved sources, after subtracting the NXB and the un- 
resolved CXB, and correcting for the telescope vignetting 
and the exposure of each pointing. This image is useful to 
set a threshold in terms of the X-ray flux. The M FLAT was 
created in the 0.7-7 keV band since the NXB is dominant 
in the higher energy band. In order to eliminate sources 
at the rim region effectively, events within a radius of 22 



mm (~ 22') from the optical axis of each GIS sensor were 
accumulated, which is by 10% larger than the radius of 
20 mm (~ 20') when we created an energy spectrum. A 
single M FLAT image was created for each sample field, by 
summing up all pointings for both GIS 2 and GIS 3 in the 
common sky coordinates. 

The NXB image to be subtracted was generated for 
each pointing and for each sensor, as described in subsec- 
tion 2.2. With regard to the unresolved CXB image, we 
created a template image for each GIS sensor, which was 
commonly used for every field and every pointing. This 
template was made from a sunlit-Earth image in 1-3 keV 
integrated over two years from 1993 June to 1995 June, 
because of its uniformness and extreme brightness. Since 
its spectrum is much softer than the CXB, we corrected 
the radial brightness profile for the XRT vignetting utiliz- 
ing a radial profile of the superposed LSS image in 0.7-7 
keV, in which irregularities due to the discrete sources are 
sufficiently smeared out. The normalization of the tem- 
plate was determined to give 88% intensity of the LSS 
field. Although this is an empirical factor, it is considered 
to be close to the unresolved fraction of the X-ray emission 
in the LSS field. Details are described by Ishisaki (1996). 

After subtracting these two images, the observed im- 
age was cross-correlated with the PSF of the XRT+GIS 
system. Since the shape of the PSF of the XRT strongly 
depends on the source position on the detector, we calcu- 
lated it by interpolating the Cyg X-l images, which were 
observed with various offset angles and azimuth angles 
within a radius of 17' from the optical axis (Takahashi et 
al. 1995; Ikebe et al. 1997). The vignetting of the XRT 
was also corrected here. This series of processes was made 
for each sensor (2 sensors) and each pointing (Np point- 
ings), and afterwards all 2 x Np images were summed up 
to build a single mosaic image. Finally, corrected for the 
exposure time considering the overlap of multiple point- 
ings, the M FLAT image was complete. 

For the thus-obtained M FLAT image, we settled on a 
certain threshold level which corresponded to ~ 2 x 10 -13 
erg cm -2 s _1 (0.7-7 keV), and searched for peaks above 
it. This flux level was chosen because the sources were de- 
tected at a 5-10 a significance, depending on the distance 
from the optical axis, for a typical exposure of 30-40 ks. 
If we assumed a power-law spectrum of T ~ 1.7, which 
is typical for the resolved sources, the threshold flux did 
not change very much (less than 10%) by altering the en- 
ergy range from 0.7-7 keV into 2-10 keV. When a peak 
was detected, a circular region was masked out from the 
data. The radius of the mask was determined so that 
the remaining surface brightness due to the tail of the 
PSF would become less than 10% of the unresolved CXB. 
This calculation was conducted using equation (1), and 
the typical mask diameters were 9' and 14' for fluxes of 
2x 10~ 13 erg cm~ 2 s _1 and lx 10~ 12 erg cm -2 s^ 1 , respec- 
tively By omitting the masked sky regions, we created an 
energy spectrum for the remained area by accumulating 
events in each field. Since the outermost region of the 
GIS f.o.v. is dominated by the NXB, we collected events 
within a radius of 20 mm (~ 20') from the optical axis of 
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each sensor. 

There are 331 sources listed in the AMSS catalog in 
our 88 selected AMSS fields. As shown in figure 5, 90% 
of those AMSS sources were eliminated in this process at 
a flux level of 4 x 10~ 13 erg cm~ 2 s _1 in 0.7-7 keV. Some 
fraction of sources fainter than the threshold of ~ 2 x 10 -13 
erg cm~ 2 s _1 (0.7-7 keV) are also eliminated because they 
happen to exist neighboring other bright sources. For each 
field, the resultant sky area after the source elimination 
is compared with the original value in figure 6. The sky 
coverage was thus reduced by about 20% with this process. 

2-4- Instrumental Response 

An accurate instrumental (XRT+GIS) response is re- 
quired to quantify the CXB spectrum. However the XRT 
response depends significantly on the off-axis angle, mak- 
ing the overall response to the CXB much different from 
that to a point source. The superposition of point-source 
responses within the integration region is still insufficient 
due to the stray light. Moreover, the responses arc slightly 
different among the 91 sample fields, since the masking 
patterns of the source elimination are different from field 
to field. 

In order to include these conditions in the response, 
we performed ray-tracing simulations for all 91 fields us- 
ing the SimARF code (Ishisaki 1996; Honda et al. 1996; 
Shibata ct al. 2001), which calculates the photon detec- 
tion efficiency of the XRT+GIS system at each energy, 
and generates a so-called Auxiliary Response File (ARF). 
In those simulations, we assumed a flat surface brightness 
which extends far beyond the GIS f.o.v., to be concrete, 
up to 2.° 5 from the optical axis of each sensor. This 
assumption is good enough as a first approximation, at 
least for the unresolved CXB after the source elimination. 
However, it leaves some room for considering that we can 
eliminate the resolved sources in the f.o.v., while we can- 
not do so for sources outside the f.o.v., from which some 
fraction of photons appear the stray light. We deal with 
this effect in subsection 5.4. 

Though a nuisance, the stray light was reproduced suf- 
ficiently well by a ray-tracing code (Tsusaka et al. 1995), 
which we calibrated against large offset observations of 
the Crab nebula up to ~ 100', as shown in figure 1. The 
major origin of the stray light was due to the X-rays to 
come through abnormal paths, e.g. only one reflection 
by the primary mirror or via reflection off the mirror 
backside, which were fully taken into account in a ray- 
tracing simulation. The stray-light estimation is consid- 
ered to be accurate to within ±10% of the CXB intensity, 
with only a mild energy dependence. We have confirmed 
that the SimARF generates an identical ARF for a point 
source in the f.o.v. with that made by ascaarf v2.81 in 
combination with the standard calibration database files, 
namely, gis2/3_ano_on_f If _180295 . f its (telescope def- 
inition file), xrt_ea_v2_0 . f its (effective area file), and 
xrt_psf_v2_0.fits (PSF file). 

A so-called ARF filter was also applied to our SimARF 
ARF, which normalized the observed flux of the Crab neb- 
ula with XRT+GIS to the previously reported level and 
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suppressed small residual structures in the spectral fit- 
ting. With the ARF filter, the Crab spectrum can be 
expressed by an absorbed power-law model with Galactic 
absorption of Nn = 2.90 x 10 21 cm -2 , a photon index of 
r = 2.09, and 2-10 keV flux of 2.16 x 10~ 8 erg enr 2 s^ 1 
(absorption not corrected) . The ARF filter was applied by 
default for the ascaarf ARF with the "arffil" parameter; 
the details are described in Fukazawa et al. (1997). As for 
the energy redistribution matrices, the released version of 
gis2/3v4_0 . rmf were used. Spectral fits were performed 
with XSPEC vlO.00 (Arnaud 1996). 

2.5. Systematic Errors 

In summary, the constituents of the observed GIS events 
in each field can be classified into the following four com- 
ponents: (1) the unresolved CXB originating in the f.o.v., 
(2) X-rays from outside of the f.o.v., (3) the resolved 
sources to be eliminated, and (4) the NXB. As an ex- 
ample, we plotted each spectrum for the whole LSS field 
(including a-LSS) in figure 7. In this figure, the spectra 
of (d) — (e), (e), (f), and (b) correspond to those classifi- 
cations, respectively. The LSS field is suitable for this 
purpose because of its large sky coverage and the un- 
biased field selection. The brightest source in the LSS 
is AX J131822+3347 with 2-10 keV flux of 1.3 x 10~ 12 
erg cm -2 s _1 . The count rate of each component in the 
0.7-7 (2-10) keV band is 0.055 (0.033), 0.052 (0.022), 
0.011 (0.005), and 0.045 (0.047) count s" 1 sensor" 1 , re- 
spectively. Since there is a comparable fraction of the 
stray light or the NXB to the unresolved CXB, accuracy 
of their estimation, as well as the statistics, determines 
the errors of the resultant CXB parameters. 

As mentioned in subsection 2.2, the reproducibility of 
the NXB is considered to be 3%, which must be carefully 
dealt with when comparing the CXB parameters from 
field to field. On the other hand, the estimation error 
of the stray light mainly acts as if the effective area of 
the XRT is over- or under-estimated in common among 
the sample fields; hence it does not matter severely in a 
ficld-to-field comparison. As mentioned in subsection 2.4, 
because the stray-light estimation is considered to be ac- 
curate to within ±10% of the CXB intensity, one should 
only pay attention to comparing our results on the abso- 
lute CXB intensity with those of other satellites. 

3. Spectral Results 

3.1. Individual Fields 

In the first place, all the 91 spectra were fitted with a 
single power-law model. The interstellar absorption JVh 
were fixed to the Galactic values, which range from 0.8 x 
10 20 cm" 2 to 22.3 x 10 20 cm" 2 (Dickey, Lockman 1990), 
and are listed in table 1. When we tried to fit the entire 
energy band (0.7-10 keV) with the free photon index, T, 
there remained a significant excess below ~ 2 keV in most 
of the 91 fields. The reduced \ 2 values showed a progres- 
sive improvement as the low-energy region was gradually 
cut out from the spectral fit. Two examples of the spec- 
tra are shown in the top panels of figure 8. Limiting the 
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fitted energy range to be 2-10 keV, an average power-law 
photon index of L hard = 1.400 ±0.008 (1 a error) with a 
standard deviation of 0.054 was obtained for the 91 sam- 
ple fields. This r hard agrees well with the previous results 
from ASCA (Gendreau et al. 1995; Miyaji et al. 1998) 
and from non-imaging measurements, such as HE AO 1 
A-2 (Marshall et al. 1980). 

The existence of a soft excess has also been noticed in 
previous spectral studies, e.g. ROSAT (Hasinger 1992), 
BBXRT (Jahoda et al. 1992), and ASCA (Gendreau et 
al. 1995). The present analysis shows that the excess 
component is clearly seen in most of the 91 fields in the 
energy range below 2 keV. Such a wide-spread detection 
of spectral data and the intensity variation of the soft 
component, as shown in figure 8, is a newly found fea- 
ture; the results would supplement the soft X-ray informa- 
tion obtained with the ROSAT all- sky survey. The varia- 
tion amplitude of the soft component is roughly 1 x 10 -8 
erg cm -2 s _1 sr _1 . This is a factor-of-thrce higher than 
the level which can be accounted for by the change of the 
Galactic Ah between the maximum and the minimum. 
To fit the 0.7-10 keV spectra in a consistent way, we need 
an additional spectral model for the soft component. For 
this, we have tried (1) power-law and (2) MEKAL mod- 
els. In both cases, the soft and the hard components were 
absorbed by the same Galactic Ah- 

For the double power-law fit, the data were unable to si- 
multaneously constrain the 2 spectral slopes in some fields. 
In fact, some of the fit showed a negative slope for the 
hard component, which seriously conflicts with the estab- 
lished value of r hard ~ 1.4. We, therefore, fixed the slope 
of the hard component at 1.4. The fit for the 91 spectra 
then gave an average power-law photon index for the soft 
component, r soft = 5.76 ±0.04 (1 a error) with a stan- 
dard deviation of 0.33. Spectral examples arc shown in 
the middle panels of figure 8, indicating a significant im- 
provement in the fit for the IRAS 19245-7245 field. The 
slope of r soft ~ 6, however, indicates a steep rise in the 
soft band, and would imply a very strong emission in the 
ROSAT band unless some low-energy cut off is present. 
This behavior may have resulted from the power-law mod- 
cling of the soft component, and it seems to complicate 
the analysis procedure. Our analysis of the Lockman Hole 
field, which has been extensively studied from other satel- 
lites, also yielded a rather steep value of r soft = G.O^oe 
(90% confidence level) compared to T ~ 2.1 obtained from 
ROSAT PSPC (Hasinger 1992). 

As the origin of the soft excess, two possibilities arc im- 
plied from ROSAT observations (Kerp 1994). One com- 
ponent is probably associated with the Galactic halo, and 
the other with the Local Hot Bubble (LHB). In Miyaji et 
al. (1998) using both the ROSAT PSPC and the ASCA 
GIS+SIS data, these components were fitted with two 
thermal models, and each plasma temperature was derived 
to be ~ 0.14 keV and ~ 0.07 keV, respectively However, 
it is hard for ASCA GIS to detect these components sep- 
arately, due to its limited sensitivity at the low-energy 
band. We therefore applied one thermal MEKAL model 
(Liedahl et al. 1995) for these soft thermal components. 



Following the fit with the double power-law model, we 
again fixed L hard at 1.4, Ah at the Galactic values and 
the metal abundance at 1 solar. The average tempera- 
ture kT of the thermal component was obtained to be 
0.39 ±0.03 keV (1 a) with a standard deviation of 0.26 
keV. Some examples of the spectral fit are shown in the 
bottom panels of figure 8. 

We then fixed the plasma temperature, kT, at sev- 
eral values of between 0.14 and 0.7 keV. The former 
value has been indicated by many previous measure- 
ments; recent examples are Miyaji et al. (1998) men- 
tioned above, and ROSA T/PSPC±rocket CCD experi- 
ments (Mendenhall, Burrows 2001). The latter value is 
implied by BeppoSAX/LECS (Parmar et al. 1999), al- 
though they say that such a high temperature was caused 
by inadequate modeling. As can be seen in the distri- 
bution of reduced \ 2 i n figure 9, the temperature range 
for the soft thermal component is not well constrained, 
mainly due to the limitation in the low-energy sensitivity 
of GIS. In order to examine the field-to-field fluctuation, 
we used a common kT fixed at the average value, 0.4 keV. 
The results of the spectral fits are listed in table 2. 

3.2. Integrated Spectrum 

The integrated spectrum for all 91 fields was con- 
structed in order to look into the fine spectral features. 
The total integration time amounted to 4.2 Ms and the 
spectrum is shown in figure 10, where the Galactic ab- 
sorption, Ah, was fixed at 4.0 x 10 20 cm -2 , which is the 
average value weighted by the exposure time. The fitting 
results for the integrated spectrum are listed in table 3, 
as well as those for the whole LSS field (including a-LSS), 
because the LSS is a good example of the so-called blank 
sky. The average r hard was obtained to be 1.411 ±0.007 
(90% confidence level). The % 2 value was 175 for 65 de- 
grees of freedom. 

The residuals suggest some systematic feature above 
8 keV, characterized by an intensity drop of < 30%. 
However, the systematic error of the NXB intensity is 
~ 3%, as mentioned earlier, and the panels in figure 10 
show that the residuals are very much reduced when the 
NXB level is varied by —3%. The cutoff feature is, there- 
fore, not significant considering the error of the NXB level. 
The night-earth spectrum in figure 7 shows a Cu-K line 
around 8 keV, which is from the gas support grids of GIS. 
However, this line is clearly subtracted in the integrated 
spectra. The cdgc-likc feature around 4.7 keV corresponds 
to the L-edge energy of xenon, which is the detector gas 
of GIS. In summary, we should say conservatively that 
the present GIS data do not indicate any significant devi- 
ation of the CXB spectrum from the nominal power-law 
spectrum in the energy range above 2 keV. 

In order to check the contribution of bright sources in 
the GIS f.o.v., integrated spectra without source elimina- 
tion were fitted similarly. These results are also shown 
in table 3. The photon index, r llard , seems not to be 
changed by applying the source elimination, while i^jard 
is significantly affected by ~ 13 % for the sample fields. 
This fraction is ~ 9% for the LSS, smaller than the for- 
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mer. This is because the AMSS fields contains brighter 
sources, which are often the target itself, than the LSS. 
Since the intensities after the source elimination agree well 
with each other, we can safely state that the source elim- 
ination worked well. 

3.3. Correlations between Spectral Parameters 

The correlations between the spectral parameters ob- 
tained in subsection 3.1 are shown in figures 11 and 12. 
There seems to be some correlation between and 
f hard w ifti a i ar g C relative amplitude in the variation of 
Fjp ft (1 cr = 52±| %). F™ {t does not show a negative cor- 
relation with iVu, which has been observed in very soft 
X-rays (e.g. Bowyer et al. 1968; Tanaka, Bleeker 1977). 
On the other hand, a clear anti-correlation is seen between 
r hard and F hard xhe spect r a l fit for individual field in- 
dicates a systematic correlation between the parameters, 
as shown in the right panel of figure 12, in the sense that 
phard b ecoraes smaller for a higher F-£ ard . However, the 
total range of the scatter is larger than the statistical er- 
ror. We, therefore, conclude there is a systematic ten- 
dency between fig 8 ™ and r hard . This correlation suggests 
that sources contributing to the fluctuation of the CXB at 
a level of 10 -13 erg cm~ 2 s _1 have relatively hard spectra. 

4. Large-Scale Distribution 

The present data are useful for looking at the all-sky 
distribution of the CXB with good sensitivity after point- 
source elimination. We examine here the large-scale dis- 
tribution of both the soft and hard components. 

4-.1. Distribution in Our Galaxy 

We first consider the spatial distribution of the soft and 
hard components in our Galaxy. The spectrum in each 
field was sorted with the Galactic latitude, b, and the 
longitude, I, and summed up in each cell. The spatial 
distributions are shown in figure 13. In this plot, the 
data were binned into 10° steps in b and 30° steps in 
Z, respectively. In plotting the distribution along Z, an 
additional condition of |6| < 60° was also applied. The 
-ZVh values were averaged over the respective angular cells. 
The instrumental responses were also summed within each 
cell. 

The soft component shows a nearly symmetric distribu- 
tion around a peak at the Galactic Center. This feature 
clearly indicates that the soft component is a Galactic 
emission, and is likely to extend further into the low- 
energy range. To examine this, we compared the soft- 
component intensity with the RASS intensity, as shown 
in figure 14. The correlation with the ROSAT- 3/4 keV 
+ 1.5 keV map (Snowden et al. 1995, 1997) is generally 
good, suggesting that these enhancements are associated 
with the Galactic bulge and the North Polar Spur. We also 
looked into some specific fields characterized by strong soft 
components. For example, in the IRAS 19254—7245 field, 
the RASS data also indicate a relatively strong flux. This 
supports the above view that the soft component has its 
origin within our Galaxy. 



In order to compare our results with previous studies, 
we fitted the observed (l,b) profiles of i^ ard and F£ ott 
with a finite radius disk model, which successfully mod- 
eled the 2-60 keV and 2-18 keV distributions observed 
with HE AO 1 A-2 (Iwan et al. 1982) or Ariel V SSI 
(Warwick et al. 1980), respectively. In this model, there 
is no Galactic X-ray emission outside of a disk of radius 
Rd, and the emission within R^ has an exponential de- 
pendence on the vertical distance above the plane charac- 
terized by a scale height parameter, h. The total X-ray 
intensity, J t ot(Z,6), is given by 



I to t(l,b)=I + E 



h/R s 
sin |6| 



1 — cxp 



tan |6| 
~h/R~ g 



in which x = cos I + y (R c \/R g ) 2 — sin 2 1, where Iq is the 
average isotropic extragalactic emission assumed in this 
model, E is the normalization constant for the Galactic 
emission, and R g is the distance to the Galactic Center. 
When we fit the observed distribution with this model, 
additional errors of 5.5% and 4.0% times i?* ard were added 
to fhard an( j F^p h ^ respectively, in order to adjust the 
reduced \ 2 unity. This procedure is justified because the 
observed extragalactic CXB intensity, itself, would scatter 
intrinsically by about this amount, as described in section 
5. 

The Ff oft distribution was fitted well by the finite ra- 
dius disk model, although there is a strong correlation be- 
tween the scale height h and the disk radius i?d, as shown 
in figure 15a. The best-fit parameters are Iq = (0.01 ± 
0.10) x 10~ 8 erg cm' 2 s~ x sr _1 , E = (1.9 ±0.9) x 10~ 8 
erg cm" 2 s" 1 sr" 1 , h = (0.19 ± 0.08) R g , and R d = (1.15 ± 
0.23) R g (errors are 1 cr), with X 2 /d.oI. = 107.3/87. Since 
Iq is consistent to be zero and £> Iq, we can say that the 
soft thermal component is almost entirely a Galactic emis- 
sion. The absorption-corrected flux of the hard power-law 
component in the 0.5-2 keV energy band is 2.07 x 10~ 8 
erg cm~ 2 s _1 sr _1 if we assume the average photon index 
of r hard = 1.412 over the whole energy range. Therefore, 
the contribution of the Galactic component in the 0.5-2 
keV band is calculated to be 49% at (I, \b\) = (0°,20°), and 
18% at \b\ = 90°, respectively. 

For the hard power-law component, we could hardly 
constrain the parameters for the i^prd distribution by fit- 
ting it with the finite radius disk model, as shown in fig- 
ure 15b. We therefore fixed the parameters at i?d = 2.8i? g 
and h = 0.73R g , which were the best-fit values for the 
HEAO 1 A-2 observation (Iwan et al. 1982). These val- 
ues are inside of the \ 2 < Xmin ± 1 region in figure 15b. We 
then obtained E = 3.7 ±2.3% of I = (5.61 ±0.13) x 10~ 8 
erg cm~ 2 s _1 sr _1 (errors are 1 cr), which agrees well 
with the HEAO 1 A-2 value of E = 3.14%. Hence, the 
Galactic component in the 2-10 keV band has a low sig- 
nificance. Using our best-fit values, the Galactic contri- 
bution in the 2-10 keV band was calculated to be 7.1% 
at (l,\b\) = (0°,20°), and 2.7% at |6| = 90°, although the 
errors are large. These estimations are again consistent 
with the very early results from Ariel V SSI (Warwick et 
al. 1980), 7.2% and 2.5%, respectively. We can also com- 
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pute the Galactic contribution to the CXB fluctuation to 
be 1 cr = 1.2%, based on the best-fit model. This is much 
smaller than the observed fluctuation of 1 a — 6.49^Qg? % 
(table 2). 

4-2. Dipole Fitting 

Aiming at detecting another spatial structure than the 
Galactic component, we next tried to fit the intensity dis- 
tribution with a dipole model, as a first-order approxi- 
mation of the spherical harmonics. The dipole intensity 
distribution is expressed by 



F x (0)-F x 



= = / cos 9, 

Fx 



(3) 



where is the angle between the center of each field 
and the assumed pole direction, / is the amplitude of the 
dipole, and Fx is the average of the observed flux. As for 
the pole direction, we searched for all directions over the 
sky with a step size of 5° . The hard and soft components 
were examined separately, and the errors in the fits were 
calculated in the same way as in the previous subsection. 
Figure 16 shows the distributions of the best-fit dipole 
amplitude in the Galactic coordinate. 

As expected from the previous subsection, the soft com- 
ponent (F x oft ) indicates a strong peak near the Galactic 
Center at (l,b) = (—5°,— 5°) with a relative amplitude of 
72 ±13% (1 cr error). The hard component (F x ,ard ) also 
suggests some dipole emission with a maximum ampli- 
tude of 3.0 ± 1.6% (1 a error) at (l,b) = (15°, 0°). There 
are slight shifts from the Galactic Center to the positions 
of the maximum dipole amplitude, towards the south di- 



rection for the Fv oft and east for the F- 



hard 



This is pre- 



sumably due to the influences of the brightest fields, the 
IRAS 19254 field at (l,b) = (322.°4,-28.°7) for the F x oft , 
and the PHL 5200 field at (59. °1, -49. °6) for the F£ ard . 
When these fields are removed in the dipole fittings, the 
directions of the shift change. 

Since the Galactic contribution is too much dominant, 
we further tried a dipole fitting after subtracting the best- 
fit dipole model with its pole directed to the Galactic 
Center. Then, the residuals were consistent with no dipole 
emission at the 90% confidence level in the whole sky, for 
both F x oft and F x ald . According to the Cosmic Microwave 
Background (CMB) dipole (Lineweaver et al. 1996), a 
dipole amplitude of 0.42% towards (l,b) = (264°, 48°) 
is expected due to the so-called Compton-Getting (CG) 
effect. The CG dipole originates in our local motion 
with respect to the distant X-ray frame, which is ex- 
pected to agree with the direction and velocity deter- 
mined from the CMB dipole. There is some evidence 
that the dipolar emission also exists in the X-ray band, 
e.g. Plionis and Georgantopoulos (1999) with the 1.5 keV 
map of the ROSAT all-sky survey, Scharf et al. (2000) 
with the HE AO 1 A- 2 all-sky survey in 2-10 keV. Our 
dipole amplitudes in this direction after the subtraction 
is -8.1 ±8.1% (1 cr error) for the F x oft and -0.5 ±1.0% 
(1 cr error) for the F x ard . It is almost impossible to detect 
such a small level of the dipole (/ = 0.42%) with our data, 
because the standard deviations of the observed intensity 



are much higher (52% and 6.5% for the F x oft and F-£ ard , 
respectively) and the number of sample fields are limited 
to 91. 

5. CXB Intensity and Fluctuations 

The present survey consists of 91 fields distributing 
in the whole sky, and provides an excellent database. 
Regarding the field-to-field fluctuation of the hard com- 
ponent (F x ard ), we evaluated the contribution from our 
Galaxy to be 1.2% in the previous section, which is much 
smaller than the observed fluctuation of 1 a = 6. 49+0^ % 
(table 2). Therefore, the origin of the CXB is consid- 
ered to be mostly due to the extra-galactic one, and its 
distribution reflects a beam-to-beam fluctuation of X-ray 
sources in the observed field. In this section, we will try to 
constrain the log AT- log S relation of X-ray sources in the 
flux range, S<2x 10" 13 erg cm" 2 s" 1 at 2-10 keV. We 
compare this fluctuation with simulations, and constrain 
the log A- logS 1 source distribution. 

5.1. Log A- Log S Relation 

In order to examine the observed fluctuation compared 
with the expected level from an assumed log A- log S rela- 
tion, a detailed simulation including all of the instrumen- 
tal characteristics is essential. In the ASCA data, there 
is a significant flux contribution (~ 35-45%; see subsec- 
tion 2.5 from outside of the f.o.v. in the observed CXB 
intensity. Inside the f.o.v., discrete sources brighter than 
Sw2x 10 -13 erg cm -2 s _1 (2-10 keV) were eliminated, 
whereas we could not exclude brighter sources outside of 
the f.o.v. We only picked up extremely bright sources, 
which may influence the CXB intensity by more than 2.5% 
based on the HE AO 1 A-l catalog and the RASS-BSC, 
and discarded the associated fields (subsection 2.3). We 
conducted a number of simulations to effectively constrain 
the intensity distribution and spectral properties of the in- 
cident sources, including all of these systematic effects. 

We first assumed a certain log A- log S relation which 
defines the intensity distribution of X-ray sources in the 
sky. The differential form of the flux (S) vs. number (n) 
relation is expressed using a normalization k and a slope 
7 as 



n(S) 



dN(> S) 
dS 



{S<S min ), 
mm < s), 



(4) 



where the notation N(> S) means the number density of 
sources brighter than S per steradian. The 7 value equals 
2.5 in the Euclidean Universe, which is a good approxi- 
mation, at least in the range of S > 10 -13 erg cm -2 s _1 
(2-10 keV). Then, the integrated form is given as 



f fc 0-7+1 
N(>S) = * 1 ~ 1 min 



7+1 (iSmin < S). 



(5) 



The CXB intensity in the flux range of S < So is calculated 

as 

1 0-7+2 c-7+2^ 
j _ 2 ^ min 



F(S ) = [ ° n(S)SdS = 



(6) 
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We tentatively adopt the value of k = ko = 1.58 x 10 15 , 
and 7 = 70 = 2.5, which gives a reasonable source density 
of N(> So) =N = 1.18 x 10 4 sr" 1 at S = S = 2.0 x 10" 13 
erg cm" 2 s -1 (2-10 kcV). The S m i n value is deter- 
mined to be 2.52 x 10 -15 erg cm~ 2 s _1 (2-10 keV), using 
F(S ) = F = 5.59 x 10" 8 erg cm" 2 s" 1 sr" 1 (2-10 keV) 
which matches the observed value (subsection 5.4). In 
figure 17a, the assumed log A- \ogS relation with these 
tentative parameters are shown by the solid line (Model 
I)- 

5.2. Simulation 

We next generated a number of skies, where point 
sources are randomly spread out up to the off-axis an- 
gle of 2.°5, with their source density and flux distribution 
following the assumed log AT- logS" relation. Figure 17b 
shows an example of the simulated sky. Among the simu- 
lated skies, we discarded ones with bright sources residing 
outside of the GIS f.o.v. (r > 20 mm ~ 20'), with the 
same criterion in selecting the AMSS fields as indicated 
in figure 2. In order to imitate the source elimination, we 
further removed sources inside of the GIS f.o.v. (r < 20 
mm ~ 20') if they were brighter than So = 2.0 x 10~ 13 
erg cm~ 2 s _1 (2-10 keV). In the actual source elimina- 
tion, finite regions surrounding these eliminated sources 
were also masked out, which reduced the averaged detec- 
tion area by about 20%. On the other hand, 33 out of the 
91 sample fields consisted of multi-pointing observations, 
covering a larger sky area than that with a single point- 
ing (0.40 deg 2 ). As a result, the sky area covered by one 
sample field after source elimination, i.e. the mean value 
of the "Area 2" column of table 1, is 0.44 deg 2 on average, 
which is 10% larger than the integration area for a single 
pointing without any source elimination. 

Indeed, these finite sizes of the masking around 
bright sources or overlapping multipointing implies non- 
independent data. These effects, however, reconcile with 
a non-uniform exposure, as long as the source distribution 
is random, e.g. clustering of sources in a few arcmin scale 
can be neglected. We investigated this effect in the simula- 
tion by changing the integration area, S, on the detector. 
According to equation (11), a mild S~ ^ dependence of 
the fluctuation width is expected. Five patterns of the 
integration area were tested; namely, we masked out a 
certain range of azimuth angles, like a pie cutting, with 
opening angles of 0°, 30°, 60°, 120°, and 180°. We found 
differences in the following simulation results are small 
and within statistical errors between 0° (100 % detector 
area) and 60° (83 %) cases. Hence, we refer only to the 
results with a 100% detector area hereafter. 

We lastly generated more than ten millions of photons 
for each simulated sky following its intensity map in the 
energy range of 1-12 keV, while assuming that all of the 
sources have a common power-law energy spectrum with 
r = 1.4. The number of photons was chosen so that it 
would correspond to a 100 ks exposure with the GIS. We 
conducted a full instrumental simulation of the XRT and 
the GIS for each photon, which uses basically the same 
code with the SimARF, utilizing the ray-tracing and the 



standard response matrices of gis2/3v4_0 . rmf . Photons 
detected by the GIS were collected to make an energy 
spectrum. Simulations with GIS 2 and GIS 3 were run 
separately for a sky, and the resultant two spectra were 
summed up afterwards. We performed this instrumental 
simulation for 369 simulated skies, and generated a set 
of spectra of the same number. Each spectrum was fit 
by a single power-law model in the 2-10 keV band, with a 
common response created by the SimARF in the same way 
as described in § 2.4, assuming a flat surface brightness. 
Note that the NXB was not taken into account in the 
simulation. The contribution of the NXB reproducibil- 
ity is examined separately in the following subsection, so 
that we can evaluate the pure Poisson noise effect in the 
observed fluctuation. 

5. 3. Comparison of Results 

We thus obtained 369 set of parameters, i^ ard and 
r hard , for the assumed logiV- \ogS relation. The distribu- 
tions of _F* ard and r hard for both the observation and the 
simulation are shown in figure 18. In the following com- 
parison, we must take care that each parameter has an 
individual error of fitting, i.e. statistical error. We there- 
fore calculated a weighted average, Av(y), and its 1 a 
error, <5Av(y), for a parameter y as 



JV F JV f 

Av(*)=E£/£i 



SAv(y) = 



1 ^ [y t -Av{y)Y 



(7) 



(8) 



where i, taking between 1 and Np, denotes the field ID of 
observations or simulations, and ?/,; and o~i denote the best- 
fit value and its 1 a error, respectively Intrinsic variance 
Sd 2 (y) and its 1 a error (5Sd 2 (y), after subtracting the 
statistical error, were computed as 



JVf 



sd 2 (,) = £fc^M 2 '- 



/ 



Z^T-^/E^' 0) 



SSd 2 (y) 



{y t - Av(y)}< 



JVf 



1 4^ 



JV F 1 

( 10 ) 



Equations (7)-(10) are correct as long as both the intrin- 
sic distribution of the parameter y and the distribution of 
statistical error of yi have Gaussian shapes. This assump- 
tion seems to be a good approximation for the parameters 
after source elimination, as can be seen in the left-hand 
panel of figure 18, although this is not true for those with- 
out the source elimination where the distribution of i^} ard 
have a non-symmetric shape. 

With regard to the observed i^ard a ft er ^he source elim- 
ination, Av(F£ ard ) was calculated to be (5.85 ± 0.04) x 
10 



-s 



erg cm s sr with a standard deviation of 

Sd\F^ Td ) = 0.38±g;g| x 10~ 8 erg cm" 2 s" 1 sr" 1 , which 

corresponds to a 6.49^o'gi % fluctuation. From the simu- 
lation discussed in subsection 5.2, Av(i^ ard ) was obtained 
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to be (5.77 ±0.02) x 10" 8 erg cm" 2 s" 1 sr" 1 with fluc- 
tuation width 5.06l[!;2o %. These average values are in 
good agreement between the observation and the simu- 
lation, while the observed fluctuation is larger than that 
of the simulation. This discrepancy can be explained by 
the systematic error, mainly due to a NXB reproducibility 
of ~ 3%. From table 3, the contribution of the system- 
atic error to _f* ard is estimated to be as = 3.2%. Then, 
the simulation becomes consistent with the observation at 
1 a level, if considering the systematic error. Therefore, 
we can say that the Euclidean distribution of 7 = 2.5 is 
acceptable with the present data. 

As for the observed power-law photon index, Av(T hard ) 



and JSd (r hard ) were obtained to be 1.412± 0.007 and 



-0.006' 



respe ctively, and Av(T ar ) = 1.414 ± 0.001 which is entirely consistent with a calculation by Barcons 

^(90% 



0.055 

and ySd 2 (r hard ) = 0.009±£P were calculated for the 
simulated results. The contribution of the systematic er- 
ror to r hard is estimated to be 0.025 from table 3. We 
have also examined a contribution due to the large-scale 
anisotropy using the dipole analysis, and it is estimated 
0.011 ±0.003 at most. Therefore, the observed distribu- 
tion of r hard is wider than that of the simulation, even 
though considering the systematic error and the large- 
scale anisotropy. This discrepancy is presumably due to 
the unrealistic assumption that all the sources have a com- 
mon power-law index of T = 1.4 in the simulation. We 
investigate this effect later. 

5.4. CXB Intensity 

In order to determine the correct CXB intensity in a 
well-defined flux range, some compensation for the stray 
light is required, because the ASCA data has a signifi- 
cant flux contribution from outside of the f.o.v., where 
we cannot eliminate bright point sources. In the previous 
subsection, we showed that the assumed logiV- log S rela- 
tion is in good agreement with the observed intensity and 
fluctuation of F-£ ard . We therefore calculated the correct 
CXB intensity using the assumed logiV- log S relation. 

Spectral fits for the simulated data indicate that 



brighter than So, which by chance lie outside of the f.o.v. 

By applying the compensation factor to the observed 
CXB intensity, the unresolved CXB flux from sources 
fainter than S were reduced to (5.67 ± 0.04) x 10~ 8 
erg cm~ 2 s _1 sr -1 (2-10 keV) when averaged for the 
91 sample fields and corrected for the Galactic absorp- 
tion. On the other hand, the total CXB flux integrated 
over S -> 00 was calculated to be (6.38 ± 0.04) x 10~ 8 
erg cm -2 s _1 sr -1 . Note that this value does not include 
the soft component, and that the absolute flux determined 
with ASCA has a systematic error of about 10%, when 
compared with the recent results with Chandra or XMM- 
Newton. If we convert the latter intensity into a power-law 
normalization at 1 keV, assuming a photon index of T = 
1.4, we obtain 9.66 ±0.07 photon keV -1 cm~ 2 s _1 sr" 



-1 —2 

cm s 



sr 



(5.77 ± 0.02) x 10~ 8 erg cm" 



sr 



Meanwhile, a simple integration of the assumed logiV- 
logS* relation, i.e. equation (6), gives F(So) = Fq = 
5.59 x 10 -8 erg cm" 2 s~ x sr" 1 in the flux range S < 
S = 2.0 x 10" 13 erg cm" 2 s" 1 (2-10 keV). This is 3.2% 
less than the simulation value which includes the contri- 
bution of bright discrete sources residing outside of the 
f.o.v. by chance. If we integrate this logN- log 5* curve 
to S ^ 00, F(oo) becomes 6.29 x 10 -8 erg cm~ 2 s _1 sr -1 
(2-10 keV) which is 9.0% larger than the simulated value. 
Therefore, the compensation factor is —3.2% in the flux 
range S < Sq and +9.0% for the total CXB flux, respec- 
tively. These compensation factors do not change, even 
though the logA^- \ogS relation has a break in the flux 
range of S < S , as long as the log TV- logS* has the same 
shape at S > So, where we think the Euclidean slope of 
7 = 2.5 has been almost established. This is because the 
sources affecting to the compensation factor have fluxes 



et al. (2000) of 10.0i^ photon keV 
confidence errors) using the ASCA and BeppoSAX results. 
Though this representation is conventionally used in much 
CXB literature, this value is very susceptible to the exis- 
tence of the soft component and the power-law index of 
the hard component. For example, changing the power- 
law index, T, by only 0.01 causes a shift of normalization 
at 1 keV by as much as 16%. 

5.5. Constraint on the Log N- Log S Relation 

We then tried to constrain the acceptable range of the 
logA^- logS* relation based on the analytic dependence of 
the fluctuation width for a given set of (^,7). The S m - ln 
value was determined for each (£,7) pair, using equation 
(6) with F(So) = 5.67 x 10~ 8 erg cur 2 s" 1 sr" 1 (2-10 
keV). According to Condon (1974), the fluctuation width 
of is given analytically as 



Up 



Av 



S, 



3-7 



3-7 

min 



s 



7+2 

min 



(11) 



where fl c s represents the effective beam size of the 
XRT+GIS system. If we take the values derived from 
the simulation, i.e. ap = 5.06%, etc., il c g is calculated to 
be 0.516 dcg 2 . Using this fl c ff value and given (fc,7) pair, 
we calculated equation (11) and searched for an accept- 
able range in the (k,j) plain. In the parameter search, we 
defined x 2 (°~f) as 

2 



a , J Sd 2 - (er F • Av) z - (a s • Av)^ 

X (0~F) 



<5Sd 2 



- to • Avr 



(12) 



where as denotes the systematic error of 3.2%. 

Since the % 2 value is determined by 7 and fc, namely, the 
degree of freedom equals two, the 90% confidence range 
of or can be evaluated as x 2 (°'f) < 4.61. In figure 19a, we 
plot the region where the (fc,7) pair gives the 90% confi- 
dence range of ctf, as well as that of 68.3%, 95% and 99%. 
In this plot, we converted k into A^(> So), which represents 
the N(> S) value for given (fc, 7 ) at S = S = 2 x 10~ 13 
erg cm -2 s . The resultant 90% acceptable region of the 
logiV- log S relation is shown in figure 19b by the outer 
solid lines. Roughly speaking, this region is constrained by 
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two factors. One is a constraint by the fluctuation, which 
is expressed by the thin solid lines. The other is a con- 
straint by the absolute CXB intensity shown by the thick 
solid lines, where the integrated source flux reach 100% 
of the absolute CXB intensity, when the power-law like 
log AT- logS 1 relation is extrapolated towards the fainter 
flux range. In practice, the constraint by the fluctuation 
is valid only in the brighter flux range where the Poisson 
noise, y/N, of the source count, N, in the f.o.v. is compa- 
rable to the fluctuation width, i.e., N(> S) < 50 with our 
data. On the other hand, the log AT- log S relation must 
be within the range of the constraint by the absolute CXB 
intensity, as long as the logiV- \ogS curve has a form that 
is gradually flattening. 

Although our obtained region is fairly wide, we can fur- 
ther constrain the region by rejecting an unrealistic set of 
{k, 7) which contradict with the previous results. When 
combined with the LSS and the AMSS results (Ueda et 
al. 1999a and 1999b, respectively), N(> S) value at S= S 
is determined in the range of 3.52 deg -2 < N(> So) < 
4.38 deg -2 . Recent deep observations with Chandra and 
XMM- Newton show that the flattening of the log AT- log 5 
curve comes out at a fainter flux around (1-2) xlO -14 
erg cm" 2 s" 1 (e.g. Baldi et al. 2002; Tozzi ct al. 2001). 
Hence, we made another condition of 7 < 2.5. If these two 
conditions are added, the acceptable region of the (k,j) 
pair becomes the area surrounded by a rectangle, as shown 
in figure 19a, and the acceptable log AT- logS region be- 
comes the area surrounded by the inner solid lines in fig- 
ure 19b. Therefore, the log AT- logS* relation would reach 
100% of the CXB at S < 3.7 x 10~ 15 erg cm~ 2 s" 1 if there 
were no flattening, and the source density at S = 10~ 16 
erg cm~ 2 s _1 is in the range of 1900-40000 deg -2 . 

As a confirmation, wc tried a simulation while assum- 
ing a bending log AT- logS curve, as shown in figures 17 
and 19b (Model 2). The model was acceptable with a 
fluctuation width of 5.48^0 20 We a -^ so confirmed the 
consistency between the analytic formula, equation (11), 
by running simulations for the boundary values of k and 
7. These results are summarized in table 4. 

5. 6. Spectral Distribution of Sources 

As described in subsection 5.3, the simulation could not 
reproduce the observed deviation of the CXB spectral in- 
dex, r, whereas both the average value and the deviation 
of F-hard were explained by the simulation fairly well. This 
is partly due to the assumption that all of the sources had 
a common spectral index of T = 1.4 in the simulation. 
Wc therefore investigated the dependence of the apparent 
CXB spectral index upon the intrinsic distribution of the 
source spectra. 

In order to investigate this effect, we introduced a 
Gaussian distribution of power-law index Ts for those 
sources which constitute the CXB. We assumed that the 
distribution of the source spectra is independent of the 
source flux. Wc chose these assumptions for simplicity, 
although the actual Universe is much more complex. For 
instance, hardening of source spectra towards fainter flux 
range has already been seen in the AMSS survey (Ueda 



et al. 1999b), and the hard sources are usually heavily ab- 
sorbed sources; hence, sources with intrinsically flat spec- 
tral indices arc rare. However, modeling all of these char- 
acteristics is beyond our scope. 

The simulation results for several Ts are summarized 
in table 5 along with the observed ones. The simulation 
suggests that the intrinsic deviation of the power-law in- 
dex, r, is more than 1.0, even when considering the devi- 
ation due to the systematic error of the NXB which is es- 
timated to be ^/sd 2 (r hard ) = 0.025 and the deviation due 
to the large-scale anisotropy of 0.011 ± 0.003. The intrin- 
sic source spectra must peak at significantly harder index 
of T s = 1.1 than the observed CXB index of r hard ~ 1.4 
if we suppose that their spectral deviation is about 1.0. 
This is because sources with soft spectra emit more pho- 
tons than do hard sources when they have the same flux; 
consequently, the averaged spectrum tends to be weighted 
towards the soft source. On the other hand, we also found 
that it was difficult to fit the simulated spectrum with a 
single power-law model for a larger deviation of Tg , since 
the integrated spectrum became a concave shape. These 
funny effects are in part caused by the simplistic model, 
which is discussed in subsection 6.4. 

6. Discussion 

The present ASCA observations have shown interest- 
ing features in the large-scale distribution of the CXB. 
Below, we will summarize the main results and address 
their implications for each subject. 

6.1. CXB Intensity and Energy Spectrum 

As shown in section 3, the integrated energy spectrum 
of the 2-10 keV CXB after the source elimination for the 
total of 4.2 Ms exposure is described by a power-law model 
with a photon index of T hmd = 1.411±0.007±0.025 (errors 
at 90% statistical and systematic). The slope is consistent 
with many previous measurements. When the nominal 
NXB is subtracted, the residual from the power-law model 
suggests a systematic deviation above 8 keV as if the ob- 
served spectrum has the cut off as shown in figure 10. 
However, this feature disappears when the NXB level is 
reduced by 3% within the systematic error. Therefore, we 
regard that the cut-off feature above 8 keV is not signif- 
icant. Similarly, because the other residual features are 
all within the NXB uncertainty, we conclude that the en- 
ergy spectrum of the CXB is consistent with the nominal 
power-law model in the energy range 2-10 keV. 

The intensity of the power-law component was calcu- 
lated to be 8.61 ±0.07 photon keV" 1 cm -2 s _1 sr _1 at 
1 keV after eliminating sources with a flux greater than 
S = 2 x 10 -13 erg cm -2 s _1 . If all the sources brighter 
than 2 x 10~ 13 erg cm~ 2 s _1 were included, the intensity 
would become 9.66±0.07 photon keV -1 cm~ 2 s _1 sr _1 at 
1 keV, using the compensation factor calculated in sub- 
section 5.4. This value is consistent with the ASCA SIS 
result based on a 250 ks exposure in 4 fields, i.e. 9.4 ±0.4 
photon keV -1 cm~ 2 s _1 sr" 1 (Gcndreauet al. 1995). Note 
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that these absolute fluxes with ASCA have systematic 
errors of ~ 10% because of the calibration uncertainty. 
Nevertheless, our result is considered to be the best esti- 
mate of the CXB intensity so far from two points of view: 
a well-calibrated and low-background instrument is used, 
and the source elimination (So ~2x 10~ 13 erg cm -2 s _1 ) 
and large solid angle (Q e ff ~ 50 deg 2 ) make the cosmic 
variance small enough, as expressed in equation (11). 
Compared with results from other missions, the present 
one is consistent with the BeppoSAX LECS value based 
on a two thermal plus a power-law model fit; 10. 4+}' 4 pho- 
ton keV" 1 cm" 2 s" 1 sr" 1 at 1 keV (Parmar et al. 1999). 
However, these values are much smaller than the ROSAT 
level of 13.4±0.3 photon keV -1 cm" 2 s" 1 sr" 1 at 1 keV 
(Hasingcr 1992), even if wc include the systematic error 
in the ASCA data. The relatively strong ROSAT flux is 
caused by the steeper spectrum with r ~ 2.1, which would 
be due to contamination of the soft component or the cal- 
ibration differences (Barcons et al. 2000). 



6.2. The Galactic Component 



hard 



We found that the (l,b) profiles of both Ff> ft and F x 
can be fitted with acceptable x 2 values by the finite radius 
disk model, which strongly suggests that a certain fraction 
of the X-ray emission has the Galactic origin. As for the 
soft thermal component (F^ lt ), all of the emission is con- 
sistent to be Galactic, with a scale height of h — 1.5 ±0.6 
kpc and a disk radius of Rd = 9.2 ± 1.8 kpc (errors are 
1 a), if we assume the distance to the Galactic Center to 
be i? g = 8 kpc. The volume cmissivity, r/(z g ), is expressed 
by 



AttE 



h 



■ cxp 



h 



(13) 



where z g is the height above the plane. We can compute 
the total luminosity of the soft component by integrating 

■n(z g ) as 



V(z s )dz g 



2-Krdr = 2(4:TrE)(nRi). (14) 



Hence, Lf ft = (1.2 ± 0.7) x 10 39 erg s" 1 (0.5-2 keV). 
Snowdcn et al. (1997) derived the total luminosity to be 
~ 2 x 10 39 erg s _1 from the R OS A T all-sky survey, which 
is consistent with our result. Note that the absolute value 
of Fj? is very sensitive to the kT value, which does vary 
between 0.14 and 0.7 keV, bringing in a systematic error 
by about an order in L^ lt . Assuming the MEKAL model 
of the temperature fcT = 0.4 keV and the metal abundance 
at one solar, the electron density in the Galactic plane is 
calculated as n c = (1.5 ±0.5) x 10" 3 cm" 3 . 

In the same way, the Galactic component in the hard 
band (Fjf rd ) is calculated as L^ ard - 8 x 10 38 erg s" 1 (2- 
10 keV), although its detection is marginal. The _F£ ard 
also showed a weak dipole feature with an amplitude of 
3.0 ± 1.6% and the peak position near the Galactic Center, 
even though the pointed fields all lie above |6| — 10°. This 
result remained essentially the same when we simply took 
the count-rate data above 3 keV instead of the i^! ard ob- 



tained from the spectral fit. The soft component can be 
approximated by a thermal model with kT ~ 0.4 keV, 
and even if we take the highest value of F x detected in 
the IRAS 19254-7245 field, the flux contribution above 3 
keV is only ~ 0.02 %. This suggests that the distribution 
of the 2-10 keV X-rays in our Galaxy really has a high 
scale-height component. Essentially the same results were 
previously reported by Warwick et al. (1980) and Iwan et 
al. (1982) with the Ariel Va,nd HEAO 1 A- 2 all-sky sur- 
veys, respectively. 

Kokubun (2001) recently examined the spatial distri- 
bution of the so-called Galactic bulge and the ridge X- 
ray emission, which has a broad enhancement around the 
Galactic Center and plane, and showed that the spectrum 
can be described by a mixture of two thermal components 
(~ 0.6 keV and ~ 3 keV, respectively) and a non-thermal 
power-law component (r ~ 1.8). Based on an analysis 
of the ASCA and RXTE data, he suggested the typical 
scale-height/lcngth of the bulge emission to be b <~ 2° and 
1^7°. A comparison with the present result implies that 
this emission is extending further away (|6| > 10°) from 
the Galactic Center with low surface brightness. Ebisawa 
et al. (2001) conducted a deep X-ray survey of a region in 
the Galactic plane with Chandra ACIS-I, and found that 
the hard X-ray emission from the Galactic ridge is truly 
diffuse, not resolved into discrete sources. However, there 
are difficulties in confining plasmas with such a energy 
density and temperature in the Galactic disk. The far 
extended hard X-ray emission in our data may indicate 
an escape of the high energy plasmas. The spatial dis- 
tribution determined by the ROSAT all-sky survey gives 
h ~ 0.95 kpc 1 and Rd ~ 5.6 kpc, which is slightly smaller 
than our results for the F^ . This is possibly because 
the ROSAT PSPC is sensitive in a softer band than the 
ASCA GIS. 



6.3. Log N~ Log S Relation 

Based on the observed fluctuation and the absolute in- 
tensity of the 2-10 keV CXB flux compared with the sim- 
ulation results, the acceptable log AT- log S relation was 
constrained as shown in figure 20. In the flux range above 
~ 10 -14 erg cm" 2 s" 1 , the relation follows the uniform 
Euclidean relation of dN(> S)/dS oc S~ 2 ' 5 , as already 
indicated from the previous observations, and lies on a 
fainter extension from the previous results from Ginga 
LAC (Butcher et al. 1997) and HEAO 1 A- 2 (Piccinotti 
et al. 1982). As figure 20 clearly shows, the present 90% 
confidence range of the log A- logS 1 curve is also consis- 
tent with the fainter-part results obtained from the recent 
Chandra and XMM-Newton deep-survey observations. We 
cannot say where the turn over exactly occurs, but it is 
close to the 10 -15 -10 -14 erg cm" 2 s" 1 level from the de- 
rived envelope of figure 20. Miyaji and Griffiths (2002) re- 
cently analyzed Chandra data for the Hubble Deep Field 
North with an area of 35.7 arcmin 2 by the fluctuation 
method. They derived a loose upper limit in the number 

1 In Snowdcn et al. (1997), the scale height parameter is used for 
the gas density, hence we have multiplied their scale height by 
0.5. 



No. ] 

density at the flux level of 2 x 10 -16 erg cm -2 s" 1 in the 
2-10 keV band to be less than 10000 per square degree. 
As figure 20 shows, the present ASCA result in the same 
energy band gives a consistent upper limit of ~ 22000 per 
square degree, in the sense that an extrapolation of the 
log TV- log S curve reaches 100% of the CXB intensity at 
this point. 

The recent Chandra results also suggest an ~ 8% dif- 
ference in the source contribution brighter than S > 4.5 x 
10" 16 erg cm" 2 s" 1 between the CDF-N and the CDF-S 
fields (Rosati et al. 2002). However, the measured inten- 
sity, i.e. (5.58 ±0.56) x 10" 8 erg cm" 2 s" 1 (CDF-N) and 
(5.15 ±0.49) x 10" 8 erg cm" 2 s" 1 (CDF-S) in 2-10 keV, 
have fairly large errors of ~ 10%, and both fields are re- 
stricted to relatively small areas (~ 0.1 deg 2 ). As pointed 
out by Barcons et al. (2000), the cosmic variance plays an 
important role in such small-area observations. We have 
shown in section 5 that ~ 5% of the intensity deviation 
can be reproduced even by a single log TV- log S relation, 
and that the fluctuation width is inversely proportional 
to the square root of the observed area, i.e. op oc l/-y/f2 c ff 
from equation (11). Considering that the GIS±XRT ef- 
fective beam-size, ri c ft, is about 0.5 deg 2 , it is not unusual 
that the intensities of the CDF-N and the CDF-S differ 
by ~ 8%. If we take our average value of the CXB in- 
tensity integrated up to S -> oo, i.e. (6.38 ±0.04) x 10~ 8 
erg cm -2 s" 1 sr _1 , the fraction of CXB resolved into dis- 
crete sources can be 80.8 ±7.7% in CDF-S and 87.5 ±8.8% 
in CDF-N, respectively, although these values do not take 
into account the uncertainty (~ 10%) of the absolute flux 
with ASCA. 

6.4- Spectral Fluctuation 

As shown in subsection 5.3, the observed distribution of 
r hard in dicates la = 0.055±g;3g|. If we try to explain this 
in terms of the Gaussian distribution of the spectral index 
of discrete sources, we need to set the intrinsic fluctuation 
width to be larger than 1.0 (1 a), and a significant fraction 
of the sources must have a harder spectrum than the CXB 
with an average of Ts ~ 1.1. This suggests a possibility 
that faint sources do not distribute along a simple log TV- 
log S 1 relation, but may consist of two or more different 
populations. However, the present spectral fluctuation 
may be partly coupled with the Galactic structure which 
affects the observed variation of r hard by - 0.011 ± 0.003 
at the maximum. Also, the systematic error on the index 
caused by the NXB subtraction is 0.025. The simulated 
fluctuation of r hard is ±0.031 for a single log TV- log S 
curve with an intrinsic Ts = 1.1 ± 1.0; therefore, the total 
fluctuation, including the NXB effect, amounts to 0.041 
and approaches the observed level. 

We also note that the Gaussian approximation for the 
distribution of the spectral index is not adequate. The 
actual situation may be that the distribution broadly con- 
sists of two components: one corresponds the population 
of unabsorbed AGNs with the peak of index distribution 
around Tg — 1.7, and the other contains heavily absorbed 
galaxies with their absorption ranging over a widely dif- 
ferent levels. In fact, several models have been proposed 



13 

along with this picture, e.g. by Madau et al. (1994), 
Comastri ct al. (1995), and Gilli et al. (2001), and they 
were successful in explaining the spectral shape of the 
CXB. However, the distributions of the intrinsic absorp- 
tion of AGN in these models are based on some assump- 
tions and/or limited observations of nearby AGNs (e.g. 
45 Seyferts of z < 0.025 in Risaliti et al. 1999). Because 
of the observed spectral fluctuation in our result should 
include all of the contribution from faint distant AGNs, 
it remains to be studied whether these models can also 
reproduce the spectral fluctuation with 1 a ~ 0.05 when 
observed with a solid angle of fl e g — 0.5 deg 2 . 

7. Conclusion 

We measured the absolute CXB intensity and its spec- 
trum from 50 square degrees based on the low background 
data with the ASCA GIS instrument. A total of 91 se- 
lected fields were studied after eliminating discrete sources 
with a threshold of ~ 2 x 10" 13 erg cm" 2 s" 1 (2-10 keV). 
The energy spectrum in the energy range of 0.7-10 keV 
is well described by a two component model comprizing a 
soft thermal emission with fcT~ 0.4 keV and a hard power- 
law model with photon index r hard = 1.412 ±0.007 ±0.025 
(1 a statistical and systematic errors). As for the extra- 
galactic power-law emission, the intensity fluctuation is 
1 cr = 6.49^o!6i wnuc the deviation due to the repro- 
ducibility of the NXB is 3.2%. We also observed a fluc- 
tuation in the photon index of 0.055j_^QQg, with a NXB 
contribution of 0.025. 

We detected an excess emission toward the Galactic 
Center in both the soft (0.7-2 keV) and hard (2-10 keV) 
energy bands, which could be fitted by a finite radius 
disk model. The emission extends well above \b\ = 10°. 
The soft thermal component is consistent to be totally 
Galactic, whereas, for the hard power-law component in 
the 2-10 keV band, the contribution form the Galaxy is 
3-7%. The soft component showed a strong correlation 
with the RASS map in the 3/4 keV and 1.5 keV bands. 

To understand the intrinsic properties of the CXB emis- 
sion, we carried out a simulation including the stray- 
light effect. We found the absolute CXB intensity to be 
(6.38±0.04±0.64) x 10" 8 erg cm" 2 s" 1 sr" 1 (1 a statisti- 
cal and systematic errors) in the 2-10 keV band. The ob- 
served fluctuation can be explained by the Poisson noise of 
the source count in the f.o.v. (Cl e s ~ 0.5 deg 2 ), even assum- 
ing a single log TV- log S relation on the whole sky. Based 
on the observed fluctuation and the absolute intensity, an 
acceptable region of the logTV- log 5 relation was derived, 
and the source density at S = 10~ 16 erg cm" 2 s" 1 is con- 
strained in the range of 1900-40000 deg" 2 . The results 
turned out to be consistent with previous measurements 
including the recent Chandra and XMM-Newton results. 
The spectral fluctuation was modeled with a Gaussian dis- 
tribution of photon index, and the simulation indicated 
that the intrinsic r.m.s. width of T is larger than 1.0 with 
a peak around V = 1.1, suggesting a large amount of hard 
sources and a large variation in the intrinsic source spec- 
tra. 
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Table 1. Field list. 



[Vol 



* Field 


tField 


Number of 


* Observation 




11/ 


life Exposure "Area 1 


**Area 2 


ID 


name 


pointing(s) 


date (10 


20 cm - 


- 2 )(deg) 


(deg) 


(ks) 


(deg 2 ) 


(deg 2 ) 


- 


LSS 


76 


1993/12/26 


1.1 


75.35 


83.21 


515.4 


7.19 


5.98 


1 


b-LSS 


21 


1994/01/05 


1.1 


73.86 


83.85 


130.2 


2.32 


1.86 


2 


c-LSS 


20 


1994/06/17 


1.1 


77.74 


82.61 


138.7 


2.26 


1.87 


3 


d-LSS 


18 


1995/01/07 


1.1 


80.38 


81.40 


127.4 


2.18 


1.58 


4 


DRACO 


3 


1993/06/04 


4.2 


102.28 


33.92 


57.0 


0.70 


0.66 


■5 


JUPITER 


1 


1993/06/06 


2.2 


286.22 


61.22 


20.4 


0.40 


0.32 


G 


NEP FIELD 


8 


1993/06/09 


4.2 


96.38 


29.90 


118.5 


0.58 


0.47 


7 


QSF 3 


4 


1993/07/11 


l.G 


250.83 


-51.99 


60.2 


0.43 


0.37 


8 


ARP 220 


1 


1993/07/26 


4.3 


36.48 


53.13 


23.1 


0.40 


0.39 


9 


3C 368 


1 


1993/09/12 


9.1 


37.57 


15.27 


19.9 


0.40 


0.36 


10 


IRAS F10214+4724 3 


1993/11/04 


1.2 


168.04 


55.07 


88.7 


0.56 


0.55 


11 


IRAS 09104 


1 


1993/11/12 


1.0 


180.98 


43.66 


34.5 


0.40 


0.27 


12 


MNVTEST 


2 


1993/11/29 


1.8 


323.86 


-59.29 


38.2 


0.47 


0.46 


13 


SN 1986J 


2 


1994/01/21 


7.5 


140.37 


-17.53 


81.3 


0.41 


0.28 


14 


AO 0235+164 


5 


1994/02/04 


9.0 


156.79 


-39.22 


42.4 


0.41 


0.29 


15 


NGC 1614 


1 


1994/02/16 


6.1 


204.53 


-34.47 


28.0 


0.40 


0.29 


1G 


NGC 1667 


1 


1994/03/06 


5.5 


204.13 


-30.21 


31.0 


0.40 


0.33 


17 


4C 41.17 


1 


1994/03/23 


10.4 


174.71 


17.39 


24.2 


0.40 


0.28 


18 


NGC 3690 


2 


1994/04/16 


0.9 


141.89 


55.38 


30.5 


0.54 


0.46 


19 


K416 


1 


1994/05/22 


1.8 


97.07 


62.49 


24.4 


0.40 


0.35 


20 


NGC 4449 


1 


1994/05/25 


1.4 


137.20 


72.38 


37.8 


0.40 


0.33 


21 


PHL 5200 


2 


1994/06/21 


5.2 


59.13 


-49.59 


92.7 


0.54 


0.37 


22 


NGC 4418 


1 


1994/06/11 


2.1 


289.87 


61.31 


26.8 


0.40 


0.34 


23 


GSGP 4 


4 


1994/06/24 


1.9 


260.39 


-88.46 


110.5 


1.37 


1.21 


24 


DI PEG 


1 


1994/06/25 


4.0 


96.11 


-43.64 


24.3 


0.40 


0.31 


25 


PG 1404+226 


1 


1994/07/13 


2.1 


21.30 


72.42 


27.6 


0.40 


0.28 


2G 


IRAS 15307+325 


1 


1994/07/22 


2.0 


51.95 


54.92 


30.1 


0.40 


0.39 


27 


Z SYSTEM 


5 


1994/08/09 


2.2 


209.43 


-65.08 


109.9 


1.67 


1.41 


28 


MG 2016+112 


2 


1994/10/24 


15.5 


53.88 


-14.04 


94.4 


0.79 


0.69 


29 


IRAS 20460+192 


1 


1994/10/27 


11.2 


64.55 


-14.86 


32.5 


0.40 


0.31 


30 


BLANK SKY 


1 


1994/11/10 


13.4 


89.02 


14.51 


23.6 


0.40 


0.35 


31 


EPSILON CMA 


1 


1994/11/16 


14.8 


239.85 


-11.22 


34.4 


0.40 


0.31 


32 


OJ 287 


3 


1994/11/18 


3.0 


206.80 


35.86 


88.5 


0.55 


0.41 


33 


RASS 1011+1736 


1 


1994/11/30 


3.2 


219.56 


52.58 


31.8 


0.40 


0.39 


34 


XY LEO 


2 


1994/12/02 


3.3 


217.77 


49.85 


19.0 


0.40 


0.31 


35 


Mrk 231 


3 


1994/12/05 


1.3 


122.03 


60.27 


117.1 


0.83 


0.73 


3G 


NGC 7320 


1 


1994/12/07 


7.9 


93.18 


-21.08 


28.6 


0.40 


0.30 


37 


CN LEO 


1 


1994/12/16 


2.9 


244.08 


56.23 


44.7 


0.40 


0.30 


38 


Mrk 273 


1 


1994/12/27 


1.1 


107.95 


59.61 


29.8 


0.40 


0.30 


39 


NGC 5135 


1 


1995/01/22 


4.6 


311.87 


32.50 


35.6 


0.40 


0.28 


40 


Q1508+5714 


2 


1995/03/02 


1.5 


93.25 


51.25 


63.4 


0.48 


0.37 


41 


GB 930704 


3 


1995/03/17 


5.2 


151.34 


23.89 


25.0 


1.18 


0.84 


42 


NGC 4125 


1 


1995/04/05 


1.8 


130.21 


51.27 


24.8 


0.40 


0.30 


43 


TXFS 1011+144 


1 


1995/05/18 


3.9 


224.35 


51.22 


28.1 


0.40 


0.37 


44 


MS 1054.5-0321 


1 


1995/05/23 


3.6 


256.54 


48.56 


48.6 


0.40 


0.31 


45 


CL 2236-04 


1 


1995/05/31 


4.0 


63.00 


-51.22 


25.0 


0.40 


0.35 


4G 


PG 1247+267 


1 


1995/06/17 


0.9 


273.35 


89.22 


26.1 


0.40 


0.32 


47 


NGC 4450 


1 


1995/06/20 


2.4 


273.70 


78.54 


25.6 


0.40 


0.33 


48 


QSO CLUSTER 


3 


1995/07/03 


1.1 


40.01 


79.08 


114.8 


1.19 


0.93 


49 


NGC 5084 


1 


1995/07/08 


8.2 


311.61 


40.51 


25.1 


0.40 


0.32 


50 


SA 68 


2 


1995/07/09 


4.2 


111.10 


-46.22 


47.2 


0.59 


0.56 


51 


Mrk 348 


1 


1995/08/04 


5.8 


122.28 


-30.80 


32.2 


0.40 


0.28 


52 


NGC 1332 


1 


1995/08/05 


2.2 


212.03 


-54.30 


44.4 


0.40 


0.36 


53 


PKS 0634-205 


1 


1995/10/09 


22.3 


229.89 


-12.29 


25.9 


0.40 


0.32 


54 


IRAS 20551-425 


1 


1995/10/19 


3.9 


358.50 


-40.77 


25.6 


0.40 


0.30 
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* Field 


tField 


Number of 


■'"Observation 


§ A H 




life Exposure 'Area 1 


**Area 2 


ID 


name 


pointing(s) 


date (10 


20 cm" 


- 2 )(deg) 


(deg) 


(ks) 


(deg 2 ) 


(deg 2 ) 


55 


HS 1946+7658 


1 


1995/10/21 


7.6 


109.12 


23.52 


29.7 


0.40 


0.40 


5G 


RGH 12 


1 


1995/10/29 


1.5 


190.96 


47.33 


34.9 


0.40 


0.34 


57 


UKMS 1 


1 


1995/11/12 


0.8 


156.69 


49.55 


41.0 


0.40 


0.32 


58 


F568-06 


2 


1995/11/16 


2.0 


217.50 


59.50 


38.8 


0.40 


0.29 


■59 


NGC 7217 


1 


1995/11/19 


10.5 


86.40 


-19.76 


62.0 


0.40 


0.27 


60 


NGC 1672 


1 


1995/11/23 


2.3 


268.82 


-38.88 


27.3 


0.40 


0.27 


61 


BJS 855 


1 


1995/11/27 


4.0 


250.20 


49.31 


32.2 


0.40 


0.35 


G2 


PG 1148+549 


1 


1995/12/07 


1.2 


140.39 


60.39 


27.4 


0.40 


0.31 


63 


IRAS 00317-2142 


1 


1995/12/11 


1.6 


86.77 


-83.16 


27.8 


0.40 


0.30 


64 


NGC 5005 


1 


1995/12/13 


1.1 


101.12 


79.18 


27.1 


0.40 


0.27 


65 


NGC 5018 


1 


1996/01/16 


7.0 


310.05 


43.07 


28.9 


0.40 


0.37 


66 


4C 38.41 


4 


1996/03/04 


1.0 


61.22 


42.31 


26.9 


0.44 


0.34 


67 


AR UMa 


1 


1996/04/27 


1.8 


167.63 


64.89 


40.0 


0.40 


0.36 


68 


PG 1114+445 


1 


1996/05/05 


1.8 


164.87 


64.43 


49.2 


0.40 


0.30 


G9 


HE 1104-1805 


1 


1996/05/31 


4.6 


270.78 


37.79 


30.8 


0.40 


0.28 


70 


RXJ1716. 6+670 


2 


1996/06/09 


3.7 


97.63 


34.05 


54.7 


0.40 


0.31 


71 


B2 1308+326 


2 


1996/06/10 


1.1 


86.16 


83.41 


28.0 


0.40 


0.28 


72 


Q2345+007 


2 


1996/06/26 


3.8 


92.15 


-58.00 


54.4 


0.40 


0.34 


73 


IRAS 00235+102 


1 


1996/07/07 


5.1 


112.92 


-51.58 


31.6 


0.40 


0.35 


74 


NGC 612 


1 


1996/07/12 


1.8 


261.19 


-77.05 


58.8 


0.40 


0.35 


75 


CL 0024+17 


1 


1996/07/21 


4.2 


114.54 


-45.20 


34.6 


0.40 


0.29 


76 


NGC 315 


1 


1996/08/05 


5.9 


124.57 


-32.39 


26.6 


0.40 


0.33 


77 


MS 0302.7+1658 


4 


1996/08/18 


10.9 


162.86 


-34.89 


98.7 


0.67 


0.47 


78 


V 471 TAURI 


1 


1996/08/26 


15.8 


172.44 


-27.83 


36.3 


0.40 


0.28 


79 


RX J1802+1804 


1 


1996/09/30 


8.5 


44.00 


18.82 


50.0 


0.40 


0.35 


80 


A548 


2 


1996/10/10 


1.9 


229.79 


-25.21 


62.4 


0.79 


0.59 


81 


IRAS 19254-7245 


1 


1996/10/16 


6.0 


322.40 


-28.70 


24.8 


0.40 


0.33 


82 


IRAS 07598+6508 


1 


1996/10/29 


4.3 


151.10 


32.14 


31.5 


0.40 


0.31 


83 


HERC-1 


2 


1996/11/02 


2.4 


75.56 


34.83 


45.0 


0.79 


0.59 


84 


IRAS 08572+3915 


1 


1996/11/03 


2.6 


183.33 


41.10 


22.7 


0.40 


0.31 


85 


A851 


3 


1996/11/16 


1.2 


171.64 


48.27 


81.7 


0.56 


0.48 


86 


NGC 7130 


1 


1996/11/18 


2.0 


9.78 


-50.32 


26.0 


0.40 


0.35 


87 


10303+7401 


1 


1996/11/26 


4.1 


134.66 


40.06 


32.2 


0.40 


0.36 


88 


4C 06.41 


2 


1996/12/13 


2.8 


241.14 


52.76 


67.8 


0.40 


0.28 


89 


M96 


1 


1996/12/16 


2.8 


234.51 


57.15 


30.8 


0.40 


0.39 


90 


IC 2560 


1 


1996/12/19 


6.5 


269.48 


19.12 


27.2 


0.40 


0.37 


91 


PG 0043+039 


2 


1996/12/21 


3.3 


121.04 


-58.05 


51.1 


0.80 


0.75 



* Field identification number defined in this paper. 

t Corresponding to the target name in the AMSS (Ueda et al. 2001) except for the LSS fields (Ueda et al. 1998). 
t For multi-pointing fields, the earliest observation date are listed. 

§ Galactic hydrogen column density toward to the observed field estimated from Dickey and Lockman (1990). 
II Center of f.o.v. in the Galactic coordinates. 
" Area before the source elimination. 
** Area after the source elimination. 
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Table 2. 'Results of the spectral fits. 



Field ID 


Field name 






pharcl 






^ 77'SOtt 

+x 


x 2 


d.o. 


1 


b-LSS 


1 


.1 


1 


416 +U.U32 
■^ 1D -0.032 


5.69-;- 





40 +U.09 

■ w -0.09 


88.0 


65 


2 


c-LSS 


1 


.1 


1 


4 q fi +0.031 
■^ Ou -0.032 


e rr+0.15 
o - oo -0.15 





07+0. 08 
°'-0.09 


65.7 


65 


3 


d-LSS 


1 


.1 


1 


44fi +0.036 
■^"-0.036 


c 71 +0.18 
°- ' 1 -Q.18 





07+0.10 

' -0.11 


65.9 


65 


4 


DRACO 


4 


,2 


1 


oqo+0.048 
■ Oao -0.048 


5 60+ 23 
O.OU_ 23 





rn+0.12 

,ou_ 12 


82.8 


65 


5 


JUPITER 


2 


.2 


1 


rnq+0.080 
■ duo -0.083 


r r 9 + 0.39 

°- oz -0.38 





09+0.22 

OZ -0.24 


76.5 


65 


6 


NEP FIELD 


4 


,2 


1 


090+0.037 
■ OZO -0.037 


r 70+0. 18 
t>-<°_0.18 





£9+0.09 
uz -0.08 


104.6 


65 


7 


QSF 3 


1 


.6 


1 


q^o+0.045 
■°^ z -0.045 


6 20+ 24 

°- zu -0.24 





oq+0.12 

oa -0.13 


66.5 


65 


8 


ARP 220 


4 


3 


1 


qco+0.074 
■ Odz -0.074 


r Q^+0.37 
5 - 84 -0.36 


1 


4 c+0.20 
/±o -0.22 


62.4 


65 


9 


3C 368 


9 


1 


1 


4f)7 +0.079 
Wl -0.078 


6 25+ - 40 


1 


4fi +0.23 
^ D -0.22 


96.6 


65 


10 


IRAS F10214+4724 


1 


,2 


1 


o/iq+0.035 
■°^ a -0.035 


6 41+ 019 

°-^ i -0.19 





r 7+O.II 

°'-0.10 


98.2 


65 


11 


IRAS 09104 


1 


.0 


1 


91 r+0.075 

■ zld -0.075 


7 00+ - 47 

' -UU_ .46 





/iq+0.19 
'^"^-O^l 


48.0 


65 


12 


MNVTEST 


1 


.8 


1 


4 qc;+0.056 
■^ oo -0.055 


r 01 +0.28 

°- oi -0.27 





fi 9+0.15 
DZ -0.17 


67.7 


65 


13 


SN 1986J 


7 


.5 


1 


c- 4n +0.044 
■°^ u -0.047 


6 46+ 25 

°-* D -0.25 





QO+0.15 


75.6 


65 


14 


AO 0235+164 


9 


,0 


1 


07Q+O.O66 
■°' a -0.069 


z. qr^+0.36 
°- ao -0.36 





-, 7+O.I6 
1( -0.15 


54.9 


65 


15 


NGC 1614 


6 


.1 


1 


qqc+0.083 
■ Oad -0.077 


5.43lg;» 





r 7 +0.18 
°'-0.20 


91.7 


65 


16 


NGC 1667 


5 


5 


1 


41 n+0.068 
■^ lu -0.069 


5.99±g:it 





qq+0.20 

aa -0.18 


90.5 


65 


17 


4C 41.17 


10 


4 


1 


c-c-O+0.084 
■ dda -0.085 


r. 4 c+0.40 
°-^°-0.38 





97+0.21 
z '-0.19 


73.0 


65 


18 


NGC 3690 





9 


1 


oqo+0. 077 
■ JJO -0.077 


c r-q+0.38 
°- oa -0.37 





/I9+0.18 
/±z -0.20 


81.7 


65 


19 


K416 


1 


,8 


1 


4f)fi +0.078 
■ WD -0.079 


4.88±8;i| 





9O+0.19 
zu -0.17 


62.6 


65 


20 


NGC 4449 


1 


4 


1 


qci +0.065 
■ OU1 -0.065 


5.67lg;» 





/i 7+0. 18 
■ 4 '-0.17 


59.1 


65 


21 


PHL 5200 


5 


.2 


1 


4 o 4 +0.037 
■*°*-0.040 


7 19+0.24 

' ' iz -0.24 





7q +o.i4 

' ' a -0.12 


74.7 


65 


22 


NGC 4418 


2 


.1 


1 


01 r+0.069 

■ Old -0.069 


6 25+ 37 

°- zo -0.36 





7O+0.18 
' ,u -0.20 


70.7 


65 


23 


GSGP 4 


1 


.9 


1 


410 +0.034 
■^ iu -0.033 


6 07+ 17 

°- u ' -0.17 





4 r+0.09 

■ 40 -o.io 


63.5 


65 


24 


DI PEG 


4 


,0 


1 


o 4 r+0.077 
■°^°-0.076 


6 19+ - 41 

°- iJ -0.40 





r 9 +0.21 

dz -0.20 


80.0 


65 


25 


PG 1404+226 


2 


.1 


1 


qri +0.086 
■ Odl -0.083 


6 24+ - 47 

D- Z "i_0.45 





oc+0.25 
Sb -0.23 


73.8 


65 


2G 


IRAS 15307+325 


2 


,0 


1 


0-1 o+0. 061 
■ olo -0.060 


6 21 +0 - 33 

D - Z1 -0.17 





cq+0.17 
JO -0.16 


47.8 


65 


27 


Z SYSTEM 


2 


.2 


1 


A0 - [ +0.036 
■^ zl -0.036 


r 77+O.I8 
°-' '-0.18 





oq+0. 09 
OJ -0.10 


51.9 


65 


28 


MG 2016+112 


15 


,5 


1 


/ioc+0.038 
■^ zu -0.038 


6 29+ - 20 

°- ZJ -0.20 


1 


1 1+0.10 
11 -0.09 


64.2 


65 


29 


IRAS 20460+192 


11 


.2 


1 


^q fi +0.080 
■^°°-0.077 


r qi+0.39 
J - al -0.38 


1 


1R +0.19 
1D -0.19 


67.3 


65 


30 


BLANK SKY 


13 


.4 


1 


/iqq+0. 076 
■^ OJ -0.074 


r. OQ+0.37 
d - ao -0.36 





7O+0.18 
' u -0.17 


58.2 


65 


31 


EPSILON CMA 


14 


.8 


1 


c-or+0.077 
■ dod -0.075 


r 9 r+0.33 

°- zo -0.32 





?4 +0.17 
'^-0.16 


86.1 


65 


32 


OJ 287 


3 


,0 


1 


+0.059 
■ W( -0.058 


A qa+0-26 

^• ao -0.25 





^7+0.13 
/il -0.14 


82.5 


65 


33 


RASS 1011+1736 


3 


.2 


1 


,,00+0.061 
■^ oo -0.061 


5 66+ 30 
o.oo_ 29 





9O+0.15 
zo -0.17 


74.1 


65 


34 


XY LEO 


3 


3 


1 


4ntS +0.103 
■ 4UO -0.105 


r 99+O.48 
O- zz _0.46 





qr+0.22 
Od -0.25 


78.6 


65 


35 


Mrk 231 


1 


3 


1 


070+O.O34 
■°'°-0.033 


5 92+ 17 

d - az -0.17 





q^+0. 09 
°^-0.10 


59.3 


65 


30 


NGC 7320 


7. 


9 


1 


qcn+0.078 
• JOU -0.077 


6 20+ 40 

°- zu -0.39 


1 


q 4 +0.20 
u ^-0.19 


79.9 


65 


37 


CN LEO 


2 


.9 


1 


4fifi +0.064 
■ TOD -0.064 


r C7+0.32 

°- ' -0.31 





oq+0.18 
oa -0.17 


61.0 


65 


38 


Mrk 273 


1 


,1 


1 


QQfi+0.085 
■ Jau -0.080 


4 7P+ - 36 
'°-0.35 





71+0.19 
' ' 1 -0.21 


76.1 


65 


39 


NGC 5135 


4 


6 


1 


4CO+0.073 
•^ dz -0.072 


r C7+0.35 

°- ' -0.34 





79+O.I8 
'' z -0.20 


74.3 


65 


40 


Q1508+5714 


1 


,5 


1 


oa 7+0.051 
-0.050 


6 13+ - 27 

°- io -0.27 





4n +0.13 
■^ u -0.14 


57.5 


65 


41 


GB 930704 


5 


2 


1 


q Kfi +0.079 
• ouu -0.080 


r C9 + 0.40 

°- oz -0.39 





7C+0.19 


62.9 


65 


42 


NGC 4125 


1 


.8 


1 


079+O.O74 
■°' z -0.072 


c qc+0.38 
o - ao -0.37 





4 q+o.2i 

'^ a -0.19 


48.7 


65 


43 


TXFS 1011+144 


3 


.9 


1 


ocq+0.069 
• JUJ -0.068 


r ? q + 0.35 
°- ' a -0.22 





oq+0.17 

oa -0.18 


70.6 


65 


44 


MS 1054.5-0321 


3 


6 


1 


009+O.O55 
• OOZ -0.055 


6 47+ - 31 

-0.30 





qr+0.15 
od -0.16 


70.2 


65 


45 


CL 2236-04 


4 


.0 


1 


q/LQ+0 078 
•°^ a -0.077 


r yn + 0.39 
O- <U_ 24 





qr+0.18 
■°°-0.20 


82.4 


65 


4G 


PG 1247+267 





,9 


1 


9OA+0.077 
■ zau -0.075 


fi Qfi+0.43 
D " 3D -0.42 





r fi +0.22 

OD -0.21 


74.0 


65 


47 


NGC 4450 


2 


4 


1 


44 c+0.077 
• 440 -0.075 


6.04t°; 4 ° 





fifi +0.21 
DD -0.23 


69.0 


65 


48 


QSO CLUSTER 


1 


,1 


1 


4fi n+0.033 
■ wu -0.032 


5.97l°;} 7 





01 +0.09 
ol -0.11 


83.9 


65 


49 


NGC 5084 


8. 


2 


1 


qqn+0.086 
• OJa -0.085 


5 qo+ - 43 


1 


O7+0.21 
u ' -0.20 


69.1 


65 


50 


SA 68 


4 


2 


1 


O7Q+0.054 
°' a -0.052 


5 85+ - 27 

°- oo -0.26 





S4 +0.13 
O4 -0.14 


56.6 


65 


51 


Mrk 348 


5. 


.8 


1 


0-1 7+0.079 
011 -0.079 


r O9+0-40 
o - oz -0.26 





qi+0.17 
O1 -0.19 


78.7 


65 


52 


NGC 1332 


2. 


.2 


1 


40 o+0.057 
• WO -0.058 


5 85+ - 30 

°- oo -0.29 





4 g+0.16 
' w -0.15 


62.0 


65 



X-Ray Background: Spectrum and its Large-Scale Fluctuation 



Table 2. (Continued) 



Field ID 


Field name 






phard 


X ^hard 




§ ^soft 


X 


> 


d.o.f. 


53 


PKS 0634-205 


22. 


3 


1 


2 oq+0.081 
■ zaa -0.079 


6.38ir 40 





QQ + 0.16 

8a -0.15 


73 


,4 


65 


54 


IRAS 20551-425 


3 


.9 


1 


c- 9zl +0.081 
■ JZ -0.080 


5.49lg;» 





r-O+0.22 
da -0.23 


63 


.7 


G5 


55 


HS 1946+7658 


7 


.6 


1 


ri O+0.063 

■ Jlo -0.063 


5.9718:1? 





°^-0.18 


48 


.7 


G5 


56 


RGH 12 


1 


,5 


1 


oc-r+0.066 
■ OJJ -0.065 


5 fi1+°- 32 





4Q+017 
1»_0.16 


59 


.0 


G5 


57 


UKMS 1 





,8 


1 


■ 4U8 -0.056 


5 S3+ - 28 





oq+0.15 
00-0.17 


G4 


1 


G5 


58 


F568-06 


2 


.0 


1 


070+0. 060 
■°'°-0.064 


fi 1fi+ - 35 

°' lu -0.33 





94+0.17 

z ^-0.19 


60 


.9 


G5 


59 


NGC 7217 


10 


,5 


1 


t- 4 o+0.063 
■ d ^°-0.062 


r 9 r+0.28 





t-9+0.15 
J^_0.14 


G3 


1 


G5 


GO 


NGC 1672 


2 


.3 


1 


,oq+0.079 
• 4JS -0.081 


c 44 +0.39 
°-^-0.38 





qq+0.20 
oo -0.23 


76 


1 


G5 


61 


BJS 855 


4 


,0 


1 


4 o fi +0.061 
■^ ou -0.061 


6.09t^2 





01 +0.16 
■ J± — 0.18 


72 


,9 


G5 


G2 


PG 1148+549 


1 


.2 


1 


441 +0.071 
■^^1_0.071 


6.08i°;p 





40 +0.20 
* u -0.23 


61 


1 


G5 


63 


IRAS 00317-2142 


1 


6 


1 


,00+0.076 
■^ Oa -0.078 







o 4 +0.21 
°^-0.24 


84 


,0 


G5 


G4 


NGC 5005 


1 


.1 


1 


970+0.078 
< a -0.077 


6.49lg;g 





rr+0.21 
°°-0.22 


72 


.6 


G5 


65 


NGC 5018 


7 


,0 


1 


070+0. 064 
■°' J -0.066 


6.04lg;» 





Q4+0.17 
a ^-0.19 


G8 


,9 


G5 


GG 


4C 38.41 


1 


.0 


1 


4 qc+0.083 
■ q:UO -0.082 


5 S7+ - 43 
-0.42 





r fi +0.24 
OD -0.23 


66 


.5 


G5 


G7 


AR UMa 


1 


,8 


1 


,79+0.054 
■^'^-0.050 


6 22 +0 ' 27 

u - zz -0.28 





1 q+0.16 
lo -0.13 


G4 


,5 


G5 


68 


PG 1114+445 


1 


.8 


1 


ooq+0.059 
• Jao -0.061 







9 r+0.16 
z>J_0.17 


78 


1 


G5 


69 


HE 1104-1805 


4 


6 


1 


+0.077 
■ qx " -0.076 


c 7^+0.39 
°-'°-0.38 





44+0.20 

^-0.21 


71 


.8 


G5 


70 


RX J1716.6+670 


3 


.7 


1 


071 +0.055 
■•5 1 i_o 056 


r 7 r+0.28 
J - ' °-Q 27 





90+013 
za -0 14 


60 


1 


G5 


71 


B2 1308+326 


1 


1 


1 


40 9 +0.081 
■ wz -0.082 


r r 7 +0.41 

°-°'-0.41 





9R+0.23 
zu -0.20 


110.8 


G5 


72 


Q2345+007 


3 


,8 


1 


4 c- 4 +0.054 
■+O+-0.053 


c 4 i+0.25 





oq+0.13 
■ Jt '— 0.14 


61 


.6 


G5 


73 


IRAS 00235+102 


5 


.1 


1 


A 1 7+0-068 
■*J- < -0.064 


6 17+ - 34 

U-J- ' -0.35 





90+0.16 
zo_ lg 


GG 


.7 


G5 


74 


NGC 612 


1 


,8 


1 


,-,7+0.052 
' -0.051 


5 R3+ - 25 

O.DO_q 25 





cc+0.14 
oo -0.13 


95 


1 


G5 


75 


CL 0024+17 


4 


.2 


1 


O/IO+0.068 
■ oq:O -0 063 


6.39+0T7 





-0 20 


G5 


,9 


G5 


7G 


NGC 315 


5 


.9 


1 


4fiO+0 076 
■ wa -0.077 


5 G4+ - 38 

°-°*-0.37 





01 +0.19 
Ji -0.21 


91 


.2 


G5 


77 


MS 0302.7+1658 


10 


,9 


1 


E fifi +0.049 
■ OOD -0.049 


4 07+0.21 

^• a ' -0.20 





01 +0.10 
ol -0.11 


93 


,6 


G5 


78 


V 471 TAURI 


15 


.8 


1 


09C+0. 075 
• JZU -0.074 


fi 74+ - 43 

O- '+-0.41 





40+0-14 


87 


.8 


G5 


79 


RX J1802+1804 


8, 


.5 


1 


ooo+0. 054 
■ ooo -0.053 


°- DO -0.30 





Q4+0-14 
a *-0.15 


GO 


,9 


G5 


80 


A548 


1 


.9 


1 


tot+0.048 
■ ooo -0.047 


r tcc+0.23 
O.OO_ 22 





41+0.15 
' 4i -0.14 


57. 


.5 


G5 


81 


IRAS 19254-7245 


6 


.0 


1 


oco+0.081 
■ OJJ -0.080 


5 S3+ 39 

°- oo -0.38 


2 


,7+0.24 
^' -0.25 


Gl 


.8 


G5 


82 


IRAS 07598+6508 


4 


.3 


1 


ti cr+0.070 
idld -0.070 


r or+0.33 

°- JO -0.32 





44+0-18 

^-0.19 


G8 


.2 


G5 


83 


HERC-1 


2 


.4 


1 


9O9+0.059 
■ zaz -0.055 


6-62±8;|t 





70+0. 18 
+ 8 -0.16 


8G 


.7 


G5 


84 


IRAS 08572+3915 


2 


.6 


1 


,1 7+O.O8I 
' -0.080 


6.03l°;g 





4fi +0.24 
*°-0.22 


58 


.6 


G5 


85 


A851 


1 


.2 


1 


40 9 +0.045 
■ wz -0.045 


t- fi9 +0.23 
°- DZ -0.22 





o n +o.ii 

ou -0.13 


70 


,4 


G5 


8G 


NGC 7130 


2 


.0 


1 


4fifi +0.077 
■™°-0.077 


r 9 r+0.36 

°- zo -0.35 





, 9 +0.21 
^ z -0.20 


67 


.2 


G5 


87 


10303+7401 


4 


.1 


1 


/L«8+0.064 
■ 4b8 -0.059 


6 lfi+°- 32 

D ' iD -0.33 





1 7+O.I8 
l'-0.17 


77 


.8 


G5 


88 


4C 06.41 


2 


.8 


1 


,r 9 +0.054 
■^ Oz -0.055 


5 fi0 +0 ' 27 

O.OU_ 26 





9C+0.14 
zo -0.14 


74 


1 


G5 


89 


M96 


2 


,8 


1 


/ioc+0.062 
■ 4Sb -0.063 


5 3fi+°- 30 

O.OO_ 29 





9^+0.18 
zu -0.16 


72 


,9 


G5 


90 


IC 2560 


6 


.5 


1 


qnc,+0.066 
■ oud -0.059 


°-'^-0.37 





+0.17 
Oi -0.17 


64 


.8 


G5 


91 


PG 0043+039 


3 


3 


1 


oqo+0.049 
■ JJJ -0.049 


c qc;+0.25 
°- JO -0.25 





q R +0.13 
oD -0.14 


86 


,7 


65 




"average 






1 41 2+0.007 


5 85 +lu)4 

o - oo -0.04 


u - oi -0.03 










"standard deviation 






n55+ 005 
u.uoo_ 006 


n 00+O.O3 
U -°^-0.04 


26+ 02 

u - ZD -0.02 









* Results from MEKAL plus power-law model fitting (see section 3). Error represents 90% confidence level. 
The temperature and abundance are fixed at kT = 0.4 keV and Z = 1.0 solar, respectively. 
1 Galactic absorption (10 20 cm -2 ) for both the power-law and the MEKAL components. Same as table 1. 



+ Flux of the power-law component (10~ 8 erg cm -2 s _1 sr _1 ) in the 2-10 keV band. 
§ Flux of the MEKAL component (10 -8 erg cm -2 s _1 sr _1 ) in the 0.5-2 keV band. 
II Weighted average Av and its 1 <r error <5Av calculated with equations (7) and (8). 

1 Standard deviation V Sd 2 and its 1 a error, yj Sd 2 ± <5Sd 2 — V Sd 2 , after the subtraction of statistical error 
calculated with equations (9) and (10). 
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Table 3. Fitting results for the integrated spectrum in the 91 sample fields and the whole LSS field (including a-LSS). 



NXB 

(%) 



(10 



jp hard 

^x 



-8 



(2-10 keV) 



erg cm 



- 1 ) (io 



TT" SO ft 

^X 



(0.5-2 keV) 



X /d-o.f. 



erg cm 



Integrated spectrum with source elimination 

+3 1.439™ (1.4 fix) 5.72^4 (5-90+ - 02 
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Table 4. Simulation results on the fluctuation width compared with the analytical formula. 
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Table 5. Simulation results on r hard compared with the observation. 
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Fig. 1. 0.7-10 keV count rate of the Crab nebula observed 
with GIS at various offset angles from the optical axis of each 
GIS sensor. Events within a radius of 20 mm (~ 20') were 
accumulated, and the CXB and the NXB counts were sub- 
tracted after being corrected for the dead time. The left panel 
is shown for a wide range of offset angles of up to ~ 100 ar- 
cmin. The observed count rates for each sensor are plotted in 
circles, as well as the expectations from the simulation plotted 
in crosses. The solid line represents an interpolated function 
using analytical formulae. The right hand panel is a close-up, 
when the Crab nebula is placed inside the GIS f.o.v. 




Fig. 3. Distribution of the 91 analyzed fields, shown with 
filled circles, in the Galactic coordinate. The open circles 
are other AMSS and LSS fields. The plotted symbols do not 
indicate the actual size of the observed field. 




Fig. 2. Fluxes of cataloged bright sources within 4° from 
the center of 100 pointed fields. The filled and open cir- 
cles correspond to the sources in the RASS-BSC (Voges et 
al. 1999) and in the HEAO 1 A-l X-ray source catalog (Wood 
ct al. 1984), respectively. The dotted line is used to discrim- 
inate the HEAO 1 sources which have large position errors, 
and the solid curve is for the ROSAT sources. Sources lying 
above these lines and outside of the GIS field (hatched region) 
can affect the measured CXB flux by more than 2.5%. And 
the relevant fields were discarded from the analysis. 
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Fig. 4. Top panel: Long-term change of the observed NXB 
count rate (count s~ 1 sensor - 1 ) in the 0.7-7 keV band accu- 
mulated within 20 mm (~ 20') from the optical axis of each 
GIS sensor. Each data point was determined from every 10 
ks exposure of the night Earth between 1993 June and 2000 
May, and plotted as a function of the mean ASCA time, which 
is defined as the elapsed time in second since 1993 January 1. 
Middle panel: observed NXB count rate divided by the esti- 
mated one using the HO + H2 monitor count. The long-term 
trend is fitted by the fourth order polynomial, which is indi- 
cated by the solid line. Bottom panel: residuals of the fitting. 
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Fig. 5. Distribution of the 331 AMSS sources in our analyzed 
fields as a function of the 0.7-7 keV flux. In the top panel, the 
solid line shows all of the AMSS sources in the 91 fields, and 
the dotted line indicates the masked-out ones in the source 
elimination. The bottom panel shows the elimination ratio 
for the AMSS sources. 
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Fig. 6. Comparison of the sky coverage before and after the 
source elimination. The horizontal and vertical axes corre- 
spond to the area before (Area 1) and after (Area 2) the 
source elimination, respectively. Three lines which limit the 
data scattering are also indicated. 
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Fig. 7. Comparison of the spectral components for the LSS 
field when accumulated within a radius of 20 mm (~ 20') from 
the optical axis of each GIS sensor. The spectra of all the 
pointings and both GIS sensors are summed up. From upper 
to lower (at 8 keV), each spectrum represents: (a) all events 
obtained in the LSS field after the data screening described in 
subsection 2.2; (b) the NXB spectrum estimated in the way 
as described in subsection 2.2; (c) the LSS spectrum after the 
NXB subtraction, i.e. (a) — (b); (d) the LSS spectrum after 
the source elimination described in subsection 2.3, in which 
the reduced fraction of 85.3% for the sky coverage is corrected; 
(e) the simulated component of the stray light which comes 
from the outside of the accumulation radius of 20 mm (~ 20') 
from the respective optical axes; (f) the summed spectrum of 
the resolved sources, i.e. (c) — (d). 
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Fig. 8. Two examples of the energy spectra after the point 
source elimination taken in (a) IRAS 00235+102 field and (b) 
IRAS 19254-7245 field. The top, middle, and bottom panels 
show spectral fits with single power-law, double power-law, 
and soft thermal plus power-law models, respectively. In each 
panel, the upper diagram shows the data (cross) fitted with 
each model (solid lines), and the lower diagram indicates the 
goodness of the fit. 
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Fig. 9. Distribution of the reduced \ 2 for spectral fits with 
three different models of the soft component. The plasma 
temperatures of the soft thermal component are indicated in 
each panel. 




Fig. 10. Integrated pulse-height spectrum of the CXB for an 
exposure time of 4.2 Ms and a sky coverage of 40 deg 2 . Each 
panel shows the data fitted with a 2-component model, the 
residual in unit of cr, and the data to model ratio. The left 
panel is for the nominal NXB subtraction, and the right one 
uses 3% reduced NXB. 
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Fig. 11. The left panel shows correlations between i^ ard 
and -Fx° ft i with the distribution histogram of _phard pitted 
in the top panel. The dashed line indicates the best-fit linear 
relation. The right panel shows the correlation between _FJ oft 
and Nn. In both panels, the error bars represent the 90% 
confidence levels. 




Fig. 12. The left panel shows correlation between i^hard anc j 
r hard , and the right panel shows examples of confidence con- 
tours for several fields. The error bars represent the 90% 
errors for a single parameter (Ax 2 = 2.70). The contours cor- 
respond to Ax 2 of 2.3, 4.6, and 9.2. 
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Fig. 14. Correlation between F^° ft and the RASS count in 
the 3/4 keV + 1.5 keV band for each sample field. The RASS 
count is averaged among the 36' X 36' region centered the 
ASCA field. 
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Fig. 15. Confidence contours between the scale height 
(h/R g ) and disk radius (Rd/Rg) parameters of the finite ra- 
dius disk model for the FJ oft and F^ ard . The contour levels 

are xLn + O- 5 ' + L0 ' + 4 - 61 ( 90 %)> + 5 - 99 (95%), +9.21 (99% 
confidence range), respectively, from inner to outer. The filled 
circle in the left panel represents the X^in position. For the 
right panel, the x 2 ■ position is out of the panel. 



Fig. 13. Distribution of the spectral parameters in our 
Galaxy based on the fits of spatially sorted pulse-height spec- 
tra. Left and right panels are for the galactic longitude (I) and 
latitude (6), respectively. The four rows from top to bottom 
show r hard , i^ ard , F-S- , and number of fields, respectively. 
The error bars correspond to the 90% confidence levels. 
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Fig. 17. Panel (a) shows two assumed log TV- logS 1 curves. 
Model 1 is the nominal curve (7 = 2.5) compared with the 
observed data, and Model 2 includes the flattening in the low 
flux range (S < 2 X 10~ 14 erg cm~ 2 s _1 ), as reported recently. 
Panel (b) shows the simulated sample image (r < 2,°5) follow- 
ing the logTV- logS relation. The central circle represents the 
CIS field of view, r = 20 mm ~ 20'. 
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Fig. 16. Distributions of dipolc amplitude for FJ (left) and i^ ard (right) plotted in the Galactic coordinate. The relative 
amplitude is indicated by color bars. The contours correspond to the 90% and 95% confidence levels from outer to inner, where the 
dipole amplitude has a non-zero value. 




Fig. 18. Distribution of the 2-10 keV flux F£ ard (left pan- 
els) and the photon index r hard (right panels). The top 
row represents the observed distribution, with solid and dot- 
ted histograms showing with and without the source elim- 
ination, respectively. The bottom panels show simulation 
results. The solid and dotted lines are for a fixed energy 
spectrum (Tg = 1.4 ±0.0) and for the Gaussian distribution, 
Tg = 1.1 ± 1.0, respectively. In both cases, the source elimi- 
nation was applied. 



Fig. 19. In panel (a), the curves indicate the acceptable 
range in the N(> So) versus 7 plane. The contours corre- 
spond to Ax 2 = 1.0 (68.3%), 4.61 (90%), 5.99 (95%), and 
9.21 (99% confidence range), respectively. The best-fit point 
is shown by a cross. The square region shows the constraint 
derived from other studies: 7 < 2.5 is given by recent Chandra 
and XMM- Newton observations (e.g. Baldi et al. 2002; Tozzi 
et al. 2001), and N(> So) = 3.95 ±0.43 deg~ 2 is estimated 
from the previous AMSS study (Ueda et al. 1999b). Panel 
(b) shows the logiV— logS relation matching with the accept- 
able JV(> So) and 7 values with two models (dotted lines). 
The outer solid lines correspond to the 90% confidence limits 
in panel (a), and the inner solid lines include the additional 
constraint by the square region in (a). In both the outer 
and the inner lines, the thick ones are boundaries, where the 
logiV- logS corves crosses the 100 % CXB flux, which we 
determined in section 5. 
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Fig. 20. logJV- logS 1 relation in the energy range 2-10 keV constrained from the present study (thick solid line), compared with a 
number of previous results. The references compiled in this plot are Miyaji and Griffiths (2002), Campana et al. (2001), Giacconi 
ct al. (2001), Mushotzky et al. (2000), Hasinger et al. (2001), Ogasaka et al. (1998), Ueda et al. (1999b), Perri and Giommi (2000), 
Gendreau et al. (1998), Butcher et al. (1997), and Piccinotti et al. (1982). 



